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Direct experimental access to some of tlie most intriguing quantum phenomena is not granted 
due to the lack of precise control of the relevant parameters in their naturally intricate environment. 
Their simulation on conventional computers is impossible, since quantum behaviour arising with 
superposition states or entanglement is not efficiently translatable into the classical language. How- 
ever, one could gain deeper insight into complex quantum dynamics by experimentally simulating 
the quantum behaviour of interest in another quantum system, where the relevant parameters and 
interactions can be controlled and robust effects detected sufficiently well. We report on the progress 
in experimentally simulating quantum many-body physics with trapped ions. 
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1. INTRODUCTION 

Simulations and deeper understanding of the dynamics 
of some tens of interacting spins are already intractable 
with the most powerful classical computers. For instance, 
the generic state of 50 spin- 1/2 particles is defined by 2^*^ 
numbers and to describe its evolution a 2^° x 2^° matrix 
has to be exponentiated [1]. Recently, one of the ten 
most powerful supercomputers, JUGENE in Jiilich, was 
exploring this regime. The current record was set by sim- 
ulating a system of 42 quantum bits (qubits) , equivalent 
to 42 spin-1/2 particles [2, 3]. 

In any case, it will not help to increase its impressive 
classical calculation capabilities to simulate only slightly 
larger quantum systems. Each doubling of the com- 
putational power will just allow to add one spin/qubit 
to the system (approximately after two years, according 
to Moore's law [4]). Furthermore, simply pursuing this 
path of exponential growth in the computer's classical 
capabilities would require an exponential shrinking of its 
electronic components[5]. The structure size currently 
amounts to approximately 30 nm, a distance spanned by 
roughly 100 atoms. If gifted engineers further minia- 
turize the sizes of their structures, "currents" of a few 
electrons will "flow" on "wires" spanned by a few atoms 
only. As a consequence, quantum effects will have to be 
considered for future classical computers, leading to se- 
rious consequences. Electrons charging a capacitor, for 
example, currently realize a storage of logical informa- 
tion: a charged capacitor represents a "one", a discharged 
capacitor a "zero". What, if the few electrons, classically 
well caught within the potential of the capacitor, follow 
their natural quantum mechanical paths and simply es- 
cape through the walls by tunnelling? 

However, allowing for quantum effects in a controlled 
way might also be exploited as a feature. Richard Feyn- 
man originally proposed [6] to use a well controlled quan- 
tum system to efficiently track problems that are very 
hard to address on classical computers and named the 
device a "quantum computer" (QC). Nowadays his pro- 
posal can be seen closer to the description of a quantum 
simulator (QS) [7]. In any case, his idea has been theoret- 
ically investigated and further developed to the concept 
of a universal QC. Fulfilling a well defined set of prerequi- 
sites, known as deVincenzo's criteria [8, 9], should make 
possible to run any classical and quantum algorithm by 
a stroboscopic sequence of operations. These have to act 
on single qubits, for example, changing their state, and 
on pairs of qubits performing changes on one qubit, con- 
ditional on the state of its mate. 

Hundreds of groups worldwide work on many ap- 
proaches in different fields of atomic, molecular and solid 
state systems to realize their version of the envisioned 
QC. For a concise review see, for example, reference 10. 

However, even assuming an ideal system and perfect 
operations will require the control of the order of 10"^ log- 
ical qubits as basis for translating any algorithm or the 
quantum dynamics of a complex system into a sequence 



of stroboscopic gate operations on a potential universal 
QC [11]. Residual decoherence will cause computational 
errors and must be minimized to allow for high oper- 
ational fidelities (99.99 ... 99.9) %) [12]. At present, 
only then, the errors could be overcome by quantum error 
correction, at the price of a reasonable but still tremen- 
dous overhead of ancilla qubits, approximately another 
100 per logical qubit. In total, of the order of 10^ qubits 
are required. Even though there appear to be no funda- 
mental obstacles for enhancing the fidelities of the oper- 
ations and for scaling the size of the systems [13[, there 
is still challenging technological development ahead. The 
realization of a universal QC is not expected within the 
next decades. 

A shortcut via analogue QS has been taken into con- 
sideration [6[ to allow to gain deeper insight into the dy- 
namics of quantum systems. "Analogue" emphasizes that 
the dynamics of the system are not translated into an al- 
gorithm of gate operations on subsets of qubits. In con- 
trast, a system of quantum particles is required, where 
(1) the initial state and its dynamics can be precisely 
controlled, (2) as many relevant parameters as possible 
manipulated and (3) the readout of the important char- 
acteristics of the final state performed in an efficient way. 
If the system's evolution was governed by a Hamiltonian 
suspected to account for the quantum effects of inter- 
est, we would be enabled to experimentally investigate 
the physics of interest isolated from disturbances, close 
to Feynman's original proposal. The requirements on 
the number of quantum particles and fidelities of opera- 
tions for analogue QS are predicted to be substantially 
relaxed compared to QC [14]. However, it remains to 
be investigated which realistic assumption on different 
sources of decoherence in the particular system will lead 
to a sufficiently small (or even appropriate) impact on 
the dedicated QS. It is predicted that QS are less prone 
to decoherence, for example, in simulating robust effects 
like quantum phase transistions (QPT). Therefore, they 
do not require any precautions in contrast to QC, which 
suffer from the costly overhead due to quantum error cor- 
rection. It was even proposed to establish decoherence as 
an asset [1]. In this context, decoherence is not to be seen 
as a source of errors, like in the field of universal QC, but 
as a resource to simulate its natural counterpart. For 
example, decoherence is suspected to be responsible and 
required for enhanced efficiencies of (quantum) processes 
in biological systems at T ^ 300 K ]15, 16]. 

To discuss the different requirements for different ana- 
logue QS, we can distinguish between two categories of 
simulations. One category deals with problems where 
QS provide a simulated counterpart that allows to ex- 
perimentally address intriguing questions that are not 
directly tractable in the laboratory. Examples are highly 
relativistic effects like Hawking and Unruh radiation or 
the Zitterbewegung of a freely moving particle predicted 
by Dirac's equation (see also reference 14). The sec- 
ond category of simulations deals with objectives that 
are (probably fundamentally) not accessible with clas- 
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sical computation, for example, the complex quantum 
dynamics of spins in solid state systems, as mentioned 
above. A promising strategy is to initialize an analogue 
QS in a state that can be prepared easily in the system of 
choice according to step (1) introduced above. Evolving 
the system adiabatically by changing its parameters ac- 
cording to (2) allows to reach a new state that is hard or 
impossible to reach otherwise, for example, via a QPT. 
The aim here is not to simulate the effects including all 
disturbances and peculiarities, because the analogue QS 
would than become as complex as the system to be simu- 
lated. The aim can be to investigate, whether the simpli- 
fied model still yields the effects observable in nature and, 
thereby, to gain a concise deeper understanding of their 
relevant ingredients. However, there remains room for 
the important discussion, whether the specific dynam- 
ics emulate nature or simulate the implemented model 
(Hamiltonian) and whether the results allow to draw fur- 
ther conclusions. 

In any case, it has to be emphasized that analogue QS 
are intrinsically not universal. That is, different realiza- 
tions of a QS will allow to simulate different systems. 
Even more important, different approaches for the iden- 
tical models (Hamiltonians) might allow to cross-check 
the validity of the QSs [17]. 

There are several systems proposed to implement ana- 
logue QS, offering different advantages [14] to address the 
physics in many-body systems. One of them consists of 
neutral atoms within optical lattices ]18, 19]. Another 
promising candidate is based on trapped ions ]20-22], 
originally suggested by I. Cirac and P. Zoller in 1995 ]23] 
in the context of QC. Trapped ions already compete at 
the forefront of many fields, were ultimate accuracy and 
precision is required, like metrology (see, for example, 
reference 24). Trapped ions offer unique operational fi- 
delities, individual addressability and short- as well as 
long-range interactions due to Coulomb forces. 

Many models of both categories of QS are promising 
candidates or already addressed by trapped ions. Ex- 
amples for the first category are emanating from the 
fields of cosmology [25-27], relativistic dynamics [28-30|, 
quantum optics [31| including quantum walks as a po- 
tential tool for QSs [16, 32-34[, chemistry [35[, and bi- 
ology [36, 37[. For the second category, quantum spin 
Hamiltonians [20[, Bose-Hubbard [21[ and spin-boson 
[38[ models were proposed to describe solid-state sys- 
tems and their simulation would allow to observe and 
investigate a rich variety of QPTs [39[. A summary and 
concise description of theoretical proposals on QS of both 
categories based on trapped ions and first experimental 
results until the year 2008 can be found in [40[. 

This reports aims to describe the current status of the 
field of experimental, analogue QS addressing many-body 
physics, its challenges and possible ways to address them. 
The first proof-of-principle experiment was achieved [41] 
and extended recently [42, 43] on a few trapped ions in 
linear radio-frequency (RF) traps. The main challenge 
for QS remains to scale up towards 50 ... 100 ions or 



even beyond. A simulated system of this size would al- 
ready reach far beyond the regime accessible via classical 
computation and, even more important, already allow to 
address open scientific questions. 

The report is organized as follows: In section 2 we 
introduce the tools available for QS by briefiy summariz- 
ing the types of traps, different ions species, and different 
technical implementations of the control of the electronic 
and motional degrees of freedom. In the following sec- 
tion 3 we derive the mathematical description based on 
references 20, 44, 45. We aim at extending the existing 
formalism to be directly applicable to more dimensions 
and individual trapping conditions envisioned in arrays 
of ions. We apply this formalism to a basic building block 
of QC, a two qubit phase gate on the radial modes mea- 
sured in our group, and emphasize similarities and dif- 
ferences between the application of similar operations for 
analogue QS. This section is supplemented by a detailed 
appendix. In section 4, we interpret the interactions in 
the context of analogue QS, which should be sufficient 
for understanding the subsequent discussion of the ex- 
perimental implementations without going through the 
details of section 3. In the following section 5, we first de- 
pict the proof-of-principle experiments on a few trapped 
ions in linear RF traps. Based on the state of the art 
capabilities we present in the second part of this section 
a summary of proposals to study many-body physics in a 
variety of solid state systems. The two following sections 
are dedicated to two proposals aiming for scaling up the 
systems: In section 6, we discuss potential realizations 
of a two-dimensional array of RF surface electrode traps. 
They are conceptually similar to promising approaches 
in Penning traps [44, 46, 47]. We will also introduce 
an alternative approach based on ions in optical traps 
in section 7, thus, trying to combine the advantages of 
trapped ions and optical lattices. Finally, we conclude in 
section 8. 



2. TOOLS REQUIRED FOR EXPERIMENTAL 
QUANTUM SIMULATIONS 

In this section, we describe the requirements to imple- 
ment analogue QS based on trapped ions. Most of these 
tools have been developed over the last decades, many 
for the purpose of quantum information processing (QIP) 
with main focus on QC. 



2.1. Scientists 

Currently, more than 30 experimental groups world- 
wide focus on QIP based on trapped ions. A rapidly 
growing percentage of those is extending their objec- 
tives into the field of analogue QS, including (to the best 
of our knowledge) groups at the University of Califor- 
nia Berkeley, Duke University, Eidgenosische Technische 
Hochschule (ETH), University of Freiburg, Gorgia Insti- 
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tute of Technology, Imperial College, University of Han- 
nover, University of Innsbruck, University of Mainz, Uni- 
versity of Maryland, Massachusetts Institute of Technol- 
ogy (MIT), Max Planck Institute for Quantum Optics 
(MPQ), National Institute of Standards and Technology 
(NIST), Sandia National Laboratories (SNL), and Uni- 
versity of Siegen. 



2.2. Ion Traps and Coulomb Crystals 

Isolating and trapping of individual particles as well as 
the precise control of their motional (external) degrees of 
freedom is key for many high precision measurements. 
Several trapping concepts have been developed for and 
implemented with ions, like RF traps [48], Penning traps 
[49] and optical traps ]50]. The physics of these devices, 
for example, of RF traps and optical dipole traps, is 
closely related. Electro-magnetic multipole fields act on 
the charge or induce electric dipole moments. The re- 
sulting forces on the particles lead in time average to a 
confining pseudopotential. The two concepts are com- 
pared in figure 1. 
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Figure 1. Two concepts for trapping charged particles. Both 
concepts require electro-magnetic fields (blue sinusoidal ar- 
row), (a) In RF traps, a RF field at frequency S7kf applied to 
quadrupole electrodes (yellow circles) interacts with a charged 
atom directly. The time averaged confining pseudopotential 
allows the ion to oscillate at frequencies ojx/y approximately 
an order of magnitude smaller than ^^rf- Since O^f > ljx/y, 
the trapping field can be understood as blue-detuned with re- 
spect to the "resonance" frequency ux / y , consequently the ion 
will seek for the field minimum in the centre of the quadrupole 
(field lines indicated by black arrows). Typical depths of the 
pseudopotential are of the order of fee x 10* K. (b) In opti- 
cal traps, the optical field is typically applied via laser beams 
that provide an intensity dependent AC Stark shift of the elec- 
tronic levels of the atom or ion. The frequency J7iaser of the 
laser can be detuned blue (red) with respect to the relevant 
electronic resonance frequency ^electron and therefore enforces 
the atom/ion to seek for low (high) fields. Typical depths of 
the pseudopotential are of the order of fee x 10"^ K. 



However, there was a delay of more than a decade be- 



tween trapping charged atoms in RF fields [51, 52] and 
trapping neutral particles with optical fields [53[ . One ex- 
planation is that RF traps provide potential depths of the 
order of several eV « /cb x 10^ K, while optical traps typi- 
cally store particles up to fee x 10~'^ K only. This discrep- 
ancy is mainly due to the comparatively large Coulomb 
force that RF fields can exert on charges. The RF field 
at typical frequencies f2RF/(27r) ~ (10 . . . 100) MHz di- 
rectly acts on the massive ion. The related motional fre- 
quencies within the deep pseudopotential amount to a 
few MHz. Optical fields, in contrast, oscillate more than 
six orders of magnitude faster: too fast for the massive 
atomic core to follow. In a simplified picture, the optical 
field has to induce a dipole moment of the electron and 
the atomic core first to allow for a subsequent interac- 
tion of the dipole with the optical field. Therefore, the 
resulting optical pseudopotential for neutral atoms and 
charged ions remains close to identical [54[. 
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Figure 2. Schematic of the three-dimensional electrode geom- 
etry of a linear RF trap, (a) Cross section through the central 
quadrupole electrodes (yellow) providing the radial confine- 
ment for the ion (blue disk), (b) Side view, where segments 
(grey) are used to apply DC voltages providing a static poten- 
tial well along the Z-ajcis. Combined with the radial {X, Y) 
pseudopotential due to the RF field, a three-dimensional con- 
finement is achieved. The ion is stored in ultra-high vacuum 
and is well protected against disturbances from the environ- 
ment. However, the fairly open geometry allows to access 
the external (motional) and internal (electronic) degrees of 
freedom, for example, with focused laser beams. 

Here we focus first on ions in linear RF traps. The 
concept for the radial confinement is depicted in figures 1 
and 2. The RF field applied to two opposing electrodes 
of the quadrupole can provide a radially confining pseu- 
dopotential. Similar to a quadrupole mass filter, one can 
find voltages for given parameters (electrode geometry 
and mass/charge ratio of the ion species) that allow for 
stable confinement in two dimensions. Additional DC 
voltages add a static harmonic potential to complete the 
three-dimensional confinement that can be assumed to 
be harmonic. Dependent on the application, these DC 
voltages can be applied to electrodes realized as rings 
or needles along the axis or by a segmentation of the 
quadrupole electrodes (see figure 2b). A confined ion 
will oscillate with frequency luz/{2tt) along the the trap 
axis and with frequencies ujx/y /i'^T^) in the radial direc- 
tions. The radial oscillation is superimposed by a fast os- 
cillation at frequency Q-rf / {2tt) (so-called micromotion), 
which increases with increasing distance of the ion from 
the trap centre, such that the RF field does not vanish 
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anymore. 

Typical parameters for conventional setups are a min- 
imal ion-electrode distance h ~ (100 . . . 1000) ym allow- 
ing for RF voltages of the order of 1000 V. 
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Figure 3. Fluorescence images of laser-cooled ions in a com- 
mon confining potential of a linear RF trap (see figure 2), 
forming differently structured Coulomb crystals, (a) A single 
ion (Mg"*-). (b) A linear string of 40 ions at ujx/y ^ ^z- 
The axis of the chain coincides with the trap Z-axis, which 
is identically orientated in the rest of the images, (c) A lin- 
ear string embedding a two-dimensional zigzag structure of 
60 ions for ujx/y > ^z- (d) A three-dimensional structure 
of more than 40 ions at ljx/y ^ '^z- The enhanced signal- 
to-noise ratio in (d) is achieved by extended exposure. Struc- 
tural phase transitions can be induced between one-, two- and 
three-dimensional crystals, for example by reducing the ratio 
of radial to axial trapping frequencies. 

Different laser cooling schemes can be applied to reduce 
the total energy of motion of the ion [45] . Doppler cooling 
[55-58] of several ions already allows to enter a regime, 
where the kinetic energy (k^T ~ mK) of the ions be- 
comes significantly smaller than the energy related to the 
mutual Coulomb repulsion. Hence, the ions cannot ex- 
change their position anymore. A phase transition from 
the gaseous (liquid) plasma to a crystalline structure oc- 
curs [59, 60]. On the one hand, the resulting Coulomb 
crystals (see figure 3) provide many similarities with solid 
state crystals already partially explaining why Coulomb 
crystals appear naturally suited to simulate many-body 
physics: (1) The ions reside on individual lattice sites. (2) 
The motion of the ions (external degree of freedom) can 
be described easiest in terms of common motional modes 



with the related quanta being phonons. The phonons 
in Coulomb crystals allow to mediate long-range inter- 
actions between the ions. In a different context, the 
phonons can also be interpreted as bosonic particles, for 
example, capable of tunnelling between lattice sites sim- 
ulated by the ions (see also section 5.2). On the other 
hand, there are advantageous differences compared to 
solid state crystals: (3) Coulomb crystals typically build 
up in ultra-high vacuum ((10~^ . . . 10~^^) mbar) and are 
very well shielded against disturbances from the environ- 
ment, thus providing long coherence times. (4) Coulomb 
crystals feature lattice constants of a few micrometres 
(see figure 3), dependent on the trapping potential coun- 
teracting the mutual Coulomb repulsion. Compared to 
a solid, where distances are of the order of Angstrom 
(10"^'^ m) the density of the structure in one dimen- 
sion is reduced by five, in three dimensions by fifteen 
orders of magnitude. This allows for individual address- 
ing of the ions and for individual preparation, control 
and readout of their electronic and motional states. (5) 
The Coulomb interaction between the charged ions is not 
shielded within the crystal, as in Coulomb crystals the 
charge of all ions has the same sign in contrast to ionic 
crystals in solid state systems. 

It has to be pointed out that it is possible to deter- 
ministically achieve phase transitions between different 
structures of Coulomb crystals for large numbers of ions 
[59, 60]. When the ratio of radial to axial confinement is 
reduced or the amount of confined ions is increased, we 
observe the transition from a linear chain of ions via a 
two-dimensional zigzag structure to a three-dimensional 
structure (see figure 3b-d). 

Despite the unique conditions in Coulomb crystals in 
linear RF traps and the high fidelities of operations, cur- 
rent experimental approaches on QS (and QC) are still 
limited to a small number of ions. The approaches in- 
clude of the order of ten ions arranged in a linear chain 
[43, 61]. This is accomplished by choosing the radial con- 
finement much stronger than the axial one. The linear 
chain orientates along the weakest (Z) direction, where 
tiny oscillations of the cooled ions around the minimum 
of the pseudopotential {X and y) and thus micromotion 
still remain negligible. 

For purposes of a QC and QS, scaling to a larger num- 
ber of spins and more dimensions while keeping sufficient 
control over all required degrees of freedom remains the 
challenge of the field. Using longer linear chains confined 
in anharmonic axial potentials [62] might provide a way 
to reach a number of ions in the system that in principle 
already exceeds capabilities of a classical supercomputer. 
Another way might be the use of RF ring traps offering 
periodic boundary conditions for static Coulomb crystals 
[59; 60] and even (more-dimensional) crystalline beams of 
ions [63-65] . A microfabricated ring trap is currently de- 
veloped and fabricated at Sandia National Laboratories 
[66]. 

The two main limitations for further scaling of the 
number of ions in a common potential are, from a practi- 
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cal point of view: (1) the emergence of 37V normal modes 
for N ions plus their sum and difference frequencies that 
lead to an increasingly crowded phonon spectrum (al- 
ready for each spatial dimension separately). Individ- 
ual spectral components become difficult to identify and 
off-resonant couplings to "spectator" transitions [44] are 
hard to avoid. However, under certain conditions, QS 
are predicted to allow for coupling to all modes simulta- 
neously, see for example reference 20. (2) QS based on 
ions in large, more-dimensional Coulomb crystals suffer 
from additional challenges, for example, intrinsic micro- 
motion (due to the displacement from the minimum of 
the pseudopotential), an inhomogeneous ion spacing (due 
to space charge effects) and the coupling between modes 
of all three spatial dimensions. 

One approach for scalability might be to generate a 
spin-off of the QIP community based on their invention 
of a new concept of a surface electrode geometry for RF 
traps [67, 68[ (see figure 13). Currently, this design is 
tested with the aim to allow for networks of intercon- 
nected linear traps. This constitutes a promising pos- 
sibility to realize the multiplex architecture of memory 
and processor traps for universal QC [69]. However, for 
QS we need a miniaturized array of traps allowing for 
more-dimensional interactions, as discussed in section 6. 

It has to be emphasized that there are other concepts 
for trapping ions, a prominent one being Penning traps. 
Penning traps provide trapping potentials of similar pa- 
rameters as RF traps. A strong, static magnetic field and 
a DC electric field yield a stable confinement of large, 
rotating Coulomb crystals. Storing many cool ions in a 
Penning trap, naturally provides, for example, a large tri- 
angular lattice of ions [70-72] that is also predicted to be 
well suited for QS [73[ . Promising results are on their way 
[46, 47]. Another challenging proposal for QS involves 
trapping ions optionally simultaneously with atoms in 
optical lattices [74[ (see section 7). 



2.3. Ions 

A huge variety of different atomic ions have already 
been used for the purpose of QIP. Every ion or, more 
specific, every isotope has different properties, for exam- 
ple, regarding the level scheme or the charge-mass ratio, 
and thus can meet different requirements of a QS. How- 
ever, they all have a single valence electron leading to 
an alkali-like level scheme. Most prominent are the earth 
alkali ions Be"*" , Mg"*" , Ca"*" , Sr"*" , and Ba''' . A similar elec- 
tronic structure possess Zn"*" , Cd^ , and Hg^ , followed by 
Yb+ [75[. 

Typically, two electronic levels with sufficiently long 
coherence times are chosen as qubit or spin states 
and It), respectively. (In principle, however, the restric- 
tion to two states is not required and the use of up to 
60 states has been proposed for (neutral) holmium [76[.) 
The types of qubits can be divided into two classes: In 
optical qubits, the states are encoded in two states with 
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Figure 4. Excerpt of the level scheme of ^^Mg"*" as an ex- 
ample of a hyperfine qubit (not to scale). ^^Mg^ has a nu- 
clear spin of / = 5/2 and thus a hyperfine-split ground state 
(Si/2,^ = 3 and 81/2,^' = 2). By applying a static mag- 
netic field of few Gauss, the degeneracy of the Zeeman sub- 
levels is lifted. The Doppler cooling laser (labeled "BD") is 
(t"*" polarized and detuned red by r/2 ~ 27r x 20 MHz from 
the cycling transition Si/2,-F = 3,Mf = 3 -f-^ P3/2,^" = 
4, Mf = 4. Here, F denotes the linewidth of the P lev- 
els. The level ||) := |Si/2, F = 3, Mf = 3) and the level 
It) •= \Si/2, F = 2, Mf = 2) are chosen as qubit states or 
(simulated) spin states, respectively. The ion is optically 
pumped into |t) during cooling. The electronic state is read 
out by a variant of "BD", which is resonant on the cycling 
transition. Hence, an ion in the state \\.) will fluoresce, while 
an ion in state |t) is off-resonant by almost 50-r and will re- 
main dark. The motional states of one of the motional modes 
are indicated as "ladders" on top of the electronic states. Two 
laser beams (labelled "Raman") detuned by A from the P3/2 
level can be used to drive two-photon stimulated- Raman tran- 
sitions between |t) and |t). A flop on the first red sideband 
is exemplary indicated by the arrows from |t) |2) — >■ |t) |1). 



a dipole-forbidden transition at an optical frequency. An 
example is ^°Ca+ with |t) := IS1/2) and |t) := ID5/2). 
The lifetime of |t) is on the order of 1 s, which defines the 
upper bound for its coherence time. In hyperfine/ Zeeman 
qubits, two sublevels from the ground state manifold are 
chosen as |t) and |t)- An applied magnetic field lifts the 
degeneracy within the manifolds of electronic levels to 
allow for spectrally resolving the dedicated states. The 
states of hyperfine/ Zeeman qubits have extremely long 
lifetimes and coherence times on the order of minutes 
have been observed [77, 78]. As an example for a hy- 
perfine/Zeeman qubit, an excerpt of the level scheme of 
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25]y[g+ ig shown in figure 4. The transition frequencies in 
hyperfine/Zeeman qubits are in the microwave regime. 

2.4. Basic Operations 

The quantized oscillation of the ions in the harmon- 
ically approximated potential of the trap gives rise to 
motional states, which are typically expressed in terms 
of Fock states Independent of the choice of qubit 
we will require three different types of couplings to elec- 
tronic states and/or motional states to assemble the tool- 
box for QC and QS based on trapped ions (for details see 
section 3). 

(a) coupling of the electronic states only (ID |n) 

It) \n)) This operation can be used to implement Rabi 
flops between the electronic states and serves as a one- 
qubit gate of a potential QC. In the context of QS it 
can be interpreted as simulated magnetic field (see also 
section 3.3). 

(b) coupling of the electronic and motional states 
(11) \n) ^ It) \n') ) This operation can drive Rabi flops 
between electronic states and different motional states, 
for example on the first red (n' = n — 1) or blue sideband 
(n' — n+1) (see also section 3.3). It can be used to create 
entanglement between the electronic and motional states 
and is an important ingredient for both sideband cooling 
and the readout of the motional state (see below). 

(c) state- dependent forces (for example, |t) \n) — )■ 
It) \n + 1)) These forces lead to state-dependent displace- 
ments. They can be used for conditional interactions 
between multiple ions, which are exploited for quantum 
gates (see sections 3.4 and 3.5) or effective spin-spin in- 
teractions in the simulation of quantum spin Hamiltoni- 
ans (see section 3.6). 

Operations (a) to (c) can be realized for both classes 
of qubits in the optical regime and for hyperfine/Zeeman 
qubits additionally via microwave fields [44, 45] : 

coupling via optical fields Optical qubit states can be 
linked by a single, nearly resonant laser beam with fre- 
quency wi and wave vector ki. Due to the lifetime of the 
It) state of 1/r ~ Is, the linewidth of the laser has to 
be very narrow 1 Hz) . Operations (a) and (b) can be 
implemented directly (see section 3.3). State-dependent 
forces (c) can be provided by a bichromatic light field 
(see, for example, references 79-81 and also section 3.4). 

In hyperfine/Zeeman qubits, the single laser beam can 
be substituted by two beams with frequencies oji, i02 and 
wavevectors ki , k2 driving two-photon stimulated- Raman 
transitions. The beams are detuned by Z\ 3> -T from 
a third level, for example, a P level (compare figure 4) 
with a typical lifetime l/F of the order of few nanosec- 
onds. In the mathematical treatment, this third level 
can be adiabatically eliminated for large detunings and 
the interaction gains the form of an interaction with a 
single beam of frequency uji = \uji — ui2\ and wavevector 
ki = ki — k^- The requirement of a narrow linewidth 
holds only for the difference frequency wi, which can be 




Figure 5. Implementations of different interaction types for 
hyperfine/Zeeman qubits. a) An operation of type (a) can be 
implemented, for example, by two-photon stimulated-Raman 
transitions driven by a pair of laser beams (shown without 
motional dependence) or directly by a microwave field. These 
types of interactions can be used for single-qubit gates in QC 
and to simulate the effective magnetic field in the simulation 
of quantum spin Hamiltonians. b) State-dependent forces (see 
type (c) in the text) can be created by two beams detuned 
by approximately the frequency of a motional mode. This 
interaction is used in the geometric phase gate for the dis- 
placement pulse [104, 117[ or in the simulation of the quantum 
Ising Hamiltonian to create the effective spin-spin interaction 
[20, 41[. 



fulfilled comparatively easily: The two beams can be gen- 
erated from the same laser using acousto-optical mod- 
ulators driven by a stable microwave source, while the 
requirements on the frequency stability of the laser are 
relaxed. For operation (a), the frequency loi has to (ap- 
proximately) meet the transition frequency of the qubit 
states (see figure 5a). For operation (b), fci in addition 
must not vanish to achieve a sufficient momentum trans- 
fer to the ions, see figure 4 and section 3.3. Therefore, 
the two beams are typically orthogonal (|fci| « v^lfcil) or 
counter-propagating (|A:i| w 2|A:i|). The state-dependent 
forces (c) can be implemented by nearly resonant beams 
{uji « 0) and beam geometries as for operation (b) (see 
figure 5b, sections 3.4 to 3.6, and section 4). 

The main technical drawback of using two-photon 
stimulated-Raman transitions is decoherence due to 
spontaneous emission after off-resonantly populating the 
third level. This limitation can be mitigated by increas- 
ing the detuning A and the intensities I1/2 of the beams, 
since the interaction strength scales with /1/2/^j while 
the spontaneous emission rate scales with 1^12! ^ ■ 

coupling via microwave fields Alternatively, transitions 
between the electronic states in hyperfine/Zeeman qubits 
can be driven laser-less by microwave fields. This al- 
lows to directly realize operation (a). However, due to 
the comparatively long wavelength and the related small 
momentum transfer {\hki\ — > 0), only negligible cou- 
pling to the motional modes can be achieved directly 
and additional efforts are required to provide opera- 
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tions (b) and (c) [82, 83]: By applying a static mag- 
netic field gradient along the axis of an ion chain, the 
transition frequency between jl) and |t) becomes site- 
dependent due to position-dependent Zeeman shifts. The 
ions can be individually addressed by applying microwave 
fields with these site-dependent frequencies ui. In addi- 
tion, this causes state-dependent forces as in the Stern- 
Gerlach experiment and allows for coupling to the mo- 
tional modes. The main challenge here is to provide 
sufficiently large magnetic field gradients and to cope 
with state-dependent transition frequencies, if high op- 
erational fidelities are required (see, for example, refer- 
ence 84). 

As an alternative to the static magnetic field gradients, 
alternating magnetic fields due to microwave currents in 
electrodes of surface electrode traps (see section 6) have 
been proposed [85] and first promising results have been 
achieved [86, 87]. Due to the small height of the ion 
above the electrode in this type of traps, a sufficiently 
large AC Zeeman shift can be generated, which can be 
treated analogously to the AC Stark shift created by laser 
beams in two-photon stimulated-Raman transitions dis- 
cussed above. However, the small height above the elec- 
trodes leads to further challenges (compare section 6) and 
high microwave powers are required. 



2.5. Initialization and Readout 

initialization of motional and electronic states The ini- 
tialization into one of the qubit states, for example |^), 
can be achieved with near-unity efficiency by optical 
pumping [88]. Regarding the motional modes, the ini- 
tialization includes Doppler cooling in all three dimen- 
sions leading to a thermal state with an average phonon 
number h of typically a few to ten quanta. This pre- 
cooling is required to reach the Lamb-Dicke regime (see 
section 3.3), where subsequent resolved sideband cool- 
ing [89-91] or cooling utilizing electro-magnetically in- 
duced transparency can be applied [92, 93]. These cool- 
ing schemes lead close to the motional ground state |0) 
(n = with probability of 98% in reference 91) of the 
dedicated modes. 

readout of electronic and motional states We dis- 
tinguish the two electronic states by observing state- 
dependent laser fiuorescence. The dipole allowed transi- 
tion to an excited state starting in the state \\.) is driven 
resonantly (see transition labelled "BD" in figure 4) in 
a closed cycle completed by spontaneous emission back 
to the state H) due to selection rules. For state |t) the 
detection laser is off-resonant. The ion therefore appears 
"bright" for while it remains "dark" for |t) [94-97]. 
Typically, a few per mill of the scattered photons are de- 
tected by a photo-multiplier tube or a CCD camera. The 
fidelity of this detection scheme has been shown exper- 
imentally to exceed 99.99 % for averaged and even indi- 
vidual measurements [98, 99]. However, additional pos- 
sibilities to enhance the detection efficiency, for example 



by methods developed for QC using ancilla qubits [100- 
102], cannot be applied to analogue QS, since all ions 
participate during the simulation. For the detection of 
the motional state, it can be mapped to the electronic 
state via an operation of type (b) and derived from the 
result of the spin state detection described above [103]. 

3. THEORETICAL EXCURSION 

The following calculations (sections 3.1 to 3.3) are the 
mathematical description of the toolbox that is available 
for both QCs and QSs. Detailed discussions can be found, 
for example, in reference 44 or reference 45. We summa- 
rize important equations in the following and extend the 
mathematical description to be applicable to scaled ap- 
proaches of QSs, for example, two-dimensional arrays of 
ions in individual traps. We will continue with the de- 
scription of an implementation of the effective spin-spin 
interaction for ions appearing in quantum spin Hamilto- 
nians. In order to investigate it isolated from other in- 
teractions, we will first discuss it from the point of view 
of quantum gates [104] in sections 3.4 and 3.5. Finally, 
in section 3.6, we use all tools to derive and discuss the 
quantum Ising Hamiltonian based on reference 20 as an 
example. The mathematical descriptions will be required 
to pursue proposals described in section 5.2. 

3.1. Theoretical Basics 

In the following we consider two-level systems only. 
The Hamiltonian describing the energy of the electronic 
states of N such systems is given by: 

Tie = E ^^i'' + . (3-1) 

—const, (omitted) 

where oj-^/i denote the energies of the states \\.) and |t), 

^ (i) 

respectively, uj^^ := ujf — lo^, and the operator az the 
Pauli operator (compare section B, equation (B.l)) act- 
ing on the i-th ion. 

The ions are considered to be trapped in a common 
harmonic potential or several individual harmonic po- 
tentials. The corresponding Hamiltonian in terms of the 
normal modes of the oscillation reads: 

/ 1\ 

m=l ^ ^ 

Here, and dj^ are the annihilation and creation oper- 
ators of the TO-th mode, respectively, and cOm the cor- 
responding frequency. In the following, the constant 
terms hjJm/2 will also be omitted and the abbreviation 
"Ho := "He + Urn will be used. 

An interaction of an ion with the electric field E oi & 
laser beam is described by —fl ■ E(f, t), where fl denotes 
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the electric dipole operator for the transition \l) o |t) 
and E{r, t) the field at the site of the ion. The Hamilto- 
nian describing the interaction of the field with N ions 
becomes: 



N 



h.c. k 



(3.3) 



Here, J?!*' = —[lE^^^^jl G M is the interaction strength 

at site i, k\ the wavevector at site i, the position 
of the i-th ion, loi the frequency of the field and (p\ an 
additional phase. In the most general form, the operator 
K^*' can be expressed as a linear combination of Pauli 
operators o^^/j^/^ ^'^'^ ^'^^ identity operator i^*^ (see sec- 
tion B): 



(3.4) 



with the prefactors cxj G M. 

The position operator r^'-* in equation (3.3) is decom- 
posed into the equilibrium position x^^ and the displace- 
ment xS^^ = f^'' — Xq*^. The terms fc/ ^ • Xq*^ give rise to a 
constant phase, which we absorb into (ySj + fci • Xq 

The displacement of the ion from its equilibrium po- 
sition x^'-* is expressed in terms of the normal modes of 
motion 

3Af 

a?^*^ = ^ {bm,i<imex + bmA+Nqmey + bmA+2N<imez) , 
m— 1 

(3.5) 

where 6,„ ^ are the elements of an (orthogonal) transfor- 
mation matrix (compare section A. equation (A.8)). Ex- 
pressing the operators of the normal modes in terms 
of the creation and annihilation operators yields 

Qrn = QmO (om + aL) with QjnO \j 2MuJ ' ''^'^^ 

where Af denotes the mass of one ion. Hence, the scalar 
product appearing in the Hamiltonian yields 



fc-w.^w^ ^^W(a„ + ai„), (3.7) 

m— 1 

where the Lamb-Dicke parameter of the m-th mode and 



3iV 



i-th site has been introduced: 

+ bm.i+Nkl ■ ey + bm,i+2Nkl ■ ez 



(3.8) 



To summarize, the interaction term of the Hamiltonian 
gains the form 



c. I kW. 



(3.9) 



3.2. Interaction Picture 



The transformation into the interaction picture 

n{ = Uliiido with f/o := e-'^"*/'* (3.10) 
can be carried out for each site separately, hence 
^^W^^Wt^W^W ^ith C/(i*) := e-'^o"*/^ (3.11) 

where is the Hamiltonian corresponding to the i-th 

ion and ?ii = X^ili ^i'''- 

The operator k*^*^ (equation (3.4)) related to the elec- 
tronic states reads as follows in the interaction picture 
(compare equation (B.15)): 



^'(0 ^ gii^t4.t*i"V2^(«)g-i'^nt*i'V2 



ai H — ^ ) c* 



(3.12) 

-Hh.c. 
(3.13) 



Here, we introduced (T^-* := ctx'^ + icTy^ and ct^!^ := 



(Jx^ — \<yy^. The terms in equation (3.9) containing the 
motional operators transform as follows (compare equa- 
tion (C.IO)): 



CXP ( '-i-Vm [aniB 



(3.14) 



Hence, it yields the following expression for the com- 
plete Hamiltonian in the interaction picture: 



H'l'^ = hnf'^ <^ exp i 



■ 3N 
m— 1 



+ h.c. U'W. 



(3.15) 



At this point fast rotating terms which average out on short timescales are neglected (rotating wave approx- 
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imation, RWA). For f]^^^ <ti uj^i we distinguish between 
two cases: In the first case, wi ^ ui^^, terms contain- 
ing e^''^'^-'-* are neglected (see figure 5b for an example 
of an implementation). (If ai = Q!2 = 0, nothing will 



change and the Hamiltonian will still be exact.) In the 
second case, \uji — oj^i\ <^ uj^i, all terms but e='='('^t^-'^i)* 
are neglected (see figure 5a). 



3iV 



/(i)(RWA) 



f?!'-* exp I i 



in—l 



3N 



(aai'^'^ + as&i'^^ + h.^ 



En 

jn—1 



(0 



(ai + ^) -7+^ + h.c. (|a;i-wt4|< 

(3.17) 



(3.16) 



3.3. iTx/ffy Interaction 



The time evolution corresponding to T^^t'^^^^"^' in 
equation (3.17) is involved. The time evolution is cal- 



I 

culated for a single ion i and a single motional mode to, 
for example, in references 44 and 45. As some simplifi- 
cations (Lamb-Dicke regime, see below) are not always 
justified for experiments, we will summarize this calcula- 
tion here. 

In this case the Hamiltonian simplifies: 



/{BWA) 



-J7iexp (i [77 (ae 



+ fltc'"*) - (cji - uj^^) t + ifi]) (ai + y ) 6-+ + h.: 



(3.18) 



Writing the state vector in the basis of electronic states 
|s) and motional Fock states 

s&{l,f} » 

the Schrodinger equation yields 

.,'l^'(«^^)|.,„)c,,„(t). 



se{4-,t} " 



(3.20) 

Matrix elements of the Hamiltonian vanish for s' = s. 
We obtain for the non- vanishing matrix elements [89, 105] 
(compare section D) 

X (n'|i)(i77e-*) \n) {t \a+\ i) 

= hQr.'.n (ai + ^) il"'-"l 
X gi([("'~")'^~('^i~'^t^)]*+Vi) 

(3.21) 



I 

where D{\) 
ator and 



j,Aat A a (jenotcs the displacement oper- 



r^ic-^'/^^i"'-"! 



lyS|^(lr"l)(,^). (3.22) 



Here, l1")(x) are the associated Laguerre polynomials, 
?i< := min(n',n), and 7i> := max(n',n). Analogously, 

we obtain {i,n\'Hi''^''\ t, ^0 = (t, ^^'1^^"^^^ | Un)* . 

We define 5 := (wi — i^^i) — [n' — 71)10. For small de- 
tunings <C w and interaction strengths |J7„'.„| ^ uj 
(resolved sideband regime), we apply a RWA neglecting 
terms rotating faster than e^"'*. Equation (3.20) can 
then be solved for each subset \n' and \n,l) separately: 



Af2„ 



(« 



1 + 



Q!2 



(3.23) 

Ci.n{t) = -inn',n (^1 + ^) ^ (-i) I"' ""I e'(^*+^^) Ct,„' (i) . 

(3.24) 

The solution of the system of differential equations yield 
Rabi oscillations between the states \l, n) ^ jt, 71') (com- 
pare section E): 
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cos(X„/,„t) + f sin(X„/,„t) 1 e 



-i5t/2 



X, 



^ sin(X„,,„0e'^*/2 



-ir2„-,„(ai + a2/i)il"'-"lc''^i and X„ 



4 I n',n I ■ 



1/2 



In the Lamb-Dicke regime, rji^ia + a^Yj ^ 1; 
equation (3.22) can be expanded to first order in "q: 

^n"-!^! = f^irj^/n (first red sideband) (3.26) 

12^^°^) = Qi (carrier) (3.27) 



a 



(LDR) 
n+l,n 



np]\/n + 1 (first blue sideband) (3.28) 



Successive red sideband transitions |i) |n) — > 
It) I?! — 1) followed by dissipative repumping to 
14,) |?7 — 1) with high probability are routinely used 
for sideband cooling close to the motional ground state 
|n 0) [89-91]. 

If the Lamb-Dicke parameter becomes effectively zero, 
the motional dependence will vanish (see equations (3.26) 
and (3.28)). The only remaining transition is the carrier 
transition equation (3.27) affecting the electronic states 
only. This is the case, for example, for two-photon stim- 
ulated Raman transitions with co-propagating beams 
or for microwave driven transitions in hyperfine qubits, 
where A;i ~ 0. In systems with more than one ion the 
ions will not be motionally coupled. That is why, equa- 
tion (3.25) also holds for each site separately in such sys- 



^sin(X„,,„t)e-'^*/2 



Ct,n'(0)' 

Ci,«(0) , 



(3.25) 



tems. 



Equation (3.25) simplifies for resonant carrier transi- 
tions (6 = 0) and a pure interaction (ai = 1 and 
a2 = 0): 



(3.29) 



where 



cos(i?/2) -ie''^ sin(t?/2)\ , . 
-ic-'^sin(i?/2) cos(t9/2) i' ^ ' 



1? := 2f2n,nt and (p '■= (pi. The rotation matrix R{tt/2, if) 
describes a 7r/2-pulse and R{Tr,(p) a 7r-pulse with phase 
if. 

3.4. Effective (g) ct^ Interaction 

We will now discuss the case of equation (3.16) with 
ai = a2 = 0. Hence, we omit the superscript of 
the Hamiltonian indicating a RWA. In the Lamb-Dicke 

r\ / 9\ 1/2 



regime, r]„ 



^ 1, the Hamiltonian can 



be expanded to first order in the Lamb-Dicke parameters 

ing terr 
uji — LUm yields 



77m'' . A subsequent RWA neglecting terms rotating faster 
than c^'^^rnt ^-^^Y^ 



/(j)(LDR) ^ ;j/2We'("'^'*+'^S") 



/(i) (LDR) (RWA) 



3N 



1 + i E (a^e-"^"-* + a. 



rn— 1 



3N 



(aoi^^^ -f aa-ri'^) +h-c. (3.31) 

(3.32) 



Note that equation (3.32) breaks up into a sum over i each. 

terms that depend on only one mode tti and one site With the excursion in section F the total time evolu- 
tion operator in the interaction picture reads 
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/(LDR)(RWA) 



(t,io) = exp 



\ L*— 1 m—l 



^^p(-iEE E (aoi« +a3aW) ® (aoi^^U ^a^i-''^ 



■i— 1 j — 1 m—l 



(3.33) 



- to) cos (^{'^ - (^1^') - sin (dmit - to) - ((^f' - ^p^)) 



The interaction described by equation (3.33) can be in- 
terpreted as follows: The first exponential function has 
the form of a displacement operator D{\) = e'*'"*"^ 
which leads to a displacement of a coherent state by A in 
phase space. Due to the (e~"^'"(*^*o) — 1) proportionality 
of the exponent, the trajectory for a coherent state of 
each mode describes a circle in phase space (or a straight 
line in the limit Sm ~ 0)- The coherent state returns to 
its initial position at times Tm = 2Trl/5m with I S N, 
where the exponent vanishes. The second exponential 
can be expanded into a tr^ (g) ct^ interaction, a inter- 
action, and a global phase. The az 'Si terms give rise 
to a geometric phase, which increases in time t, and the 
(72 terms lead to a dynamic phase [106, 107]. The area 
in phase space enclosed by the trajectory is proportional 
to these phases. 

3.5. Geometric Phase Gates 

The collective interaction of multiple ions with the 
same laser(s) has been proposed for the implementation 
of quantum gates [79, 108-111]. These gates are de- 
scribed in the z-basis by equation (3.33) and have been 
first implemented in references 104 and 112. M0lmer- 
S0rensen gates can be mathematically treated analo- 
gously in a rotated basis and are described in detail in 
references 81 and 80. Implementations are reported in 
references 113-116. 

We will exemplary discuss geometric phase gates based 



on the Gz ® (Tz terms in equation (3.33) in the following. 
They offer excellent tools to investigate a pure az (8) (Jz 
interaction required for the simulation of more involved 
Hamiltonians like quantum spin Hamiltonians. 

The interaction according to the Hamiltonian can be 
implemented [41, 104, 117] by stimulated-Raman transi- 
tions driven by two beams with wavevectors ki , k2 and 
difference frequency close to a (several) motional mode(s) 
(see figure 5b) . On average the differential AC Stark shift 
between the levels ||) and |t) caused by the two beams 
can be compensated by choosing appropriate polariza- 
tions of the beams. Still, on short timescales ^ 27r/(5,„ 
the ions experience a state-dependent force that leads to 
the above displacement in the phase spaces of the corre- 
sponding modes. 

In the original implementation of the geometric phase 
gate [104] two ^Be"'' ions are used. The state-dependent 
forces amount to Fj, = —2F^. This means that the op- 
erators (aoi^'-* -f ctzo-f') have diagonal elements 1 and 
-2, which is fulfilled for ao = -1/2 and = 3/2. The 

-(1) -(2) 

effective wavevectors ki = ki = ki — k2 point along 
the axis of the linear trap and the laser beams are de- 
tuned by ^STR = 27r X 26 kHz from the stretch (STR) 
mode. The effect of the centre-of-mass (COM) mode can 
be neglected (^com ~ IOO^str)- The ions are placed 
at the same phase of the stimulated-Raman interaction 

{^'^'=^f'=^). 

For t ~ Tg = 27r/(5sTR the time evolution operator 
equation (3.33) simplifies: 



/(LDR)(RWA) 
I 



(Tg,0)«exp -27ri 



EE 



i-l)^-Wf4^^ 



"STR 



(* 



(3.34) 



where we have used 77STR '7str ^ ~^str ^^'^ 
neglected the global phase arising from the 1^*^ (g) 
i(^) terms. The sequence of the gate is similar to 
the one in figure 7, but without the second displace- 
ment pulse D2- Ideally, the initial state Itjj) = 
lU) I'^COM = 0, nsTR = 0) is rotated to 1/2(|U) + lit) + 
|t;) + |n)) IncOM = 0, nsTR = 0) by the first i?(^/2, 7r/2) 



pulse (the phase (p of the first pulse can be chosen arbi- 
trarily). The only non- vanishing contributions arise from 
the (Tz (g) az terms for the |4,t) and |ti) states, which gain 
a geometric phase 

^STRit/n = X 4^|^ai (3.35) 

"STR 

By choosing appropriate beam intensities and thus i7i. 
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these phases equal ??sTi?,it/t4 ^ —it/2. The subsequent 
R{tt/2,tt) and R{tt/2,tt/2) pulses lead to the final Bell 
state If/i) = 1/a/2(|||) + i |tt)): which is achieved experi- 
mentally with a fidelity of F = 97% [104]. 



a) 




— jvirt) 

= -266 kHz 

i 



ROC 

1330 kHz 



2.18 MHz 



OMHz 




= -65kHz 



3.78 MHz COM/ ^ 3.384 MHz 
^ ' A ^ 3.254 MHz 



65 kHz 



OMHz 



axial 



radial 



Figure 6. Comparison between parameters of geometric phase 
gates [104] with two ions using the axial motional modes and 
radial motional modes, (a) The parameters correspond to 
the gate from reference 117. The axial centre-of-mass (COM) 
and stretch (STR) mode have a large frequency difference 
{2tt X 1.6 MHz). The detuning of the Raman beams from the 
STR mode amounts to Sstr = -2tt x 266 MHz. That is why 
the main contribution to the differential geometric phase be- 
tween |ii)/ltt) and |it)/lt4') is due to a (single) loop in the 
phase space of the STR mode. However, as already suggested 
in reference 104, the detuning from the COM mode is chosen 
to be an integer multiple of the detuning from the STR mode 
(<5coM = —5 X (5str). Hence, there is no entanglement left be- 
tween the electronic and motional modes at the gate duration 
Tg — |27r/5sTR| ~ 3.75 lis. (Note that the spin-echo sequence 
is not included in Tg.) (b) The parameters correspond to a 
phase gate on two of the radial motional modes. The ra- 
dial centre-of-mass (COM) and rocking (ROC) mode have a 
comparatively small frequency difference of only 2-k x 130 kHz. 
The detunings from both modes are chosen to have same abso- 
lute values resulting in (approximately) equal contributions to 
the acquired geometric phase from both modes. The gate du- 
ration according to the original implementation would amount 
to Tg = |27r/5sTR| = 15.4 ys. As the displacement pulse is re- 
peated in the second gap of the spin-echo sequence (compare 
figure 7) to cancel dynamic phases (compare reference 112 
and see text), the duration increases by an additional factor 
of two. 

A similar implementation of the geometric phase gate 
is reported in reference 117 based on two ^^Mg"*" ions. 
The state-dependent forces amount to = —3/2F-f 
(ao = ^1/4 and as ~ 5/4). Furthermore, the detuning 
from the STR mode amounts to Sstr = —2tt x 266 kHz 
and simultaneously the detuning from the COM mode 
6cOM = —271 X 1330kHz (compare figure 6a). Hence, 
the effect of the COM mode is also exploited for the 



gate. As the detuning from the COM mode is cho- 
sen an integer multiple of the detuning from the STR 
mode ((5coM = ^S^str), the first exponential in equa- 
tion (3.33) still becomes unity for the gate duration of 
Tg = |27r/(5sTR|- (In other words, all circular trajectories 
in all phase spaces return to their initial position for Tg.) 
As a result, there is no entanglement left between the 
electronic and motional states. 

Analogous to equation (3.35), but considering (5str < 

and (5coM > for the detunings and ?7com ~ '^com 
for the Lamb-Dicke parameters of the COM mode, the 
geometric phases yield: 
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(3.36) 
(3.37) 



By adjusting the beam intensities appropriately the dif- 
ferential phase between |H) / WX) and / It^) can be 
adjusted to fulfil <?sTR4,t/t4. ~ ^COMU/tt = '^Z^- As 



STRit/ti 



has the opposite sign compared to 



COMU/tt' 



the geometric phase gate makes use of two motional 
modes simultaneously. 

However, some of the dynamic phases from the COM 
mode do not vanish: 



3COM 
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■^I '7cOM 

"com 



(3.38) 



These phases have an absolute value of 2ao/a3 of the 
geometric phase from the COM mode and lead to a small 
deviation from the ideal state at the end of the gate. 

Compared to the original implementation in refer- 
ence 104, the gate is speeded up by approximately a fac- 
tor of 10 and the fidelity F for the Bell state exceeds 
95%. (Note that the duration of the spin-echo sequence 
is not included in Tg, because its rotations could be much 
faster and empty gaps can in principle be removed.) 

The radial motional modes are interesting, because 
they are similar to the normal modes in systems of in- 
dividual traps for each ion (compare section 6), which 
are promising candidates for scalable systems in quan- 
tum simulations. To investigate the differences between 
the axial and radial modes of motion the geometric phase 
gate with ^^Mg"*" is performed on a pair of radial modes 
(see also reference 116 for a M0lmer-S0rensen gate per- 
formed on the radial modes). 

The detunings from the COM and ROC (the equiv- 
alent to the STR mode in terms of the axial motional 
modes) are chosen to have same absolute values ^com = 
—S-Roc = 2tt X 65 kHz (see figure 6). The geomet- 
ric phases acquired on each motional mode are basi- 
cally the same as in equation (3.36), where "STR" has 
to be replaced by "ROC", and equation (3.37). (How- 
ever, the signs change due to a change of the signs of 
the detunings.) The contributions to the total differen- 
tial geometric phase between |J,4,)/|tt) and |it)/lti) due 
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to the COM and ROC mode are (approximately) equal 
now. However, the dynamic phase (analogous to equa- 
tion (3.38)) arising from the COM mode can no longer 
be neglected. 

displacement pulses 



7r/2 



Td. 



D2 

Td. 



tt/2 



Tv/2 



Figure 7. Pulse scheme of the geometric phase gate. It 
consists of a spin-echo sequence (_R(7r/2, 7r/2), /j(7r,7r/2), 
R{n /2,TT /2) pulses) with a displacement pulse (labeled "Di") 
in the first gap of the spin-echo sequence for the original im- 
plementation of the phase gate [104]. For the gate on the 
radial modes of motion a second displacement pulse (labelled 
"D2") is introduced to cancel dynamic phases from Di due to 
different absolute values of the forces on It) and |4,) (compare 
reference 112). For all gates the duration of each displace- 
ment pulses is chosen to be Td ~ |27r/(5sTR/Roc| such that 
each displacement pulse leads to a closed loop in each phase 
space. Hence, the total gate duration amounts to Tg = Td for 
the original implementation and Tg = 2Td for the gate on the 
radial modes of motion. (Note that the spin-echo sequence is 
not included in Tg.) The dashed R(n /2,tt /2 + if) analysis 
pulse is added for the measurement of the gate fidelity (see 
figure 9). 

The pulse scheme of the geometric phase gate is modi- 
fied by adding a second displacement pulse in the second 
gap of the spin-echo sequence (see figure 7 and compare 
reference 112). The intensities of the beams are now ad- 
justed for a differential geometric phases due to each dis- 
placement pulse of ^coMU/n~^ROCit/ti = '^Z^- While 
the geometric phases of both displacement pulses add up 
to tt/2, the dynamic phases cancel each other, as the tt 
pulse of the spin-echo sequence interchanges o- |tt) 
(and lit) ^ Iti))- Additionally, the more symmetric 
pulse scheme enhances the robustness of the gate against 
uncompensated differential AC Stark shifts between ||) 
and It). 

The total fluorescence from the two ions as a function 
of the total duration of the displacements 2Ti} is shown in 
figure 8. The gate duration due to the smaller detunings 
and the second displacement pulse is more than a factor 
of 8 longer than for the gate in reference 117. Still, the 
fidelity exceeds 95 % (see figure 9). 



3.6. Quantum Ising Hamiltonian 

Above we have introduced az'S^^d'z interactions that are 
used in quantum gates. In the following we will present 
a slightly different approach, in which Ising spin-spin in- 
teractions are continuously induced by means of optical 
forces. 
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Figure 8. Total fiuorescence from the two ions as a func- 
tion of the total displacement duration 2Td (compare fig- 
ure 7). The detected fluorescence signal from both ions 
amounts to approximately 7 counts/20 jis for state |H) and 
close to zero for |tt)- The duration between the displacement 
pulses is chosen to be T„ — |27r/(5coM/ROcl in the experi- 
ment (compare figure 7). At Tg ~ 30.8 ys the state \ip} ~ 
(Itt) + i lii)) |"-C0M ~ 0, TiRoc = 0) is prepared. Each data 
point represents the average of 400 measurements (squares 
and triangles) and 200 measurements (circles), respectively. 
The statistical errors are on the order of the size of the sym- 
bols. The curve is based on a fit of the time evolution of 
equation (3.33) with an additional empirical exponential de- 
cay to mimic decoherence effects. The only fit parameters are 
the fiuorescence for |H) amounting to 7.2 counts/20 ps and 
the decay constant r ~ 290 ps. The gate serves as experimen- 
tal reference for the isolated interaction strength and is not 
optimized to provide the highest gate fidelity. 



The spin-spin interaction as proposed in reference 20 
and experimentally realized in the simulation of a quan- 
tum Ising Hamiltonian in reference 41 is identical to the 
interaction iii described in section 3.4. (Note that a sim- 
ilar proposal involving the same mathematics is given in 
reference 22.) However, the quantum Ising Hamiltonian 
contains an additional (simulated) magnetic field point- 
ing in x-direction. We will adapt our notation in this 
section and split the total interaction Hamiltonian into 
the following terms: Jig denotes the term that gener- 
ates the spin-spin interaction and 'Hm denotes the term 
leading to the simulated magnetic field. The index "I" of 
the frequencies !7i and uji etc. is changed to "S" or "M" 
accordingly in the respective terms. The complete inter- 
action is described by the Hamiltonian "Hi = Us + 'Hm- 
In the following, we will first derive the spin-spin inter- 
action Hamiltonian from T-ls focusing on an Ising inter- 
action [ffz 'Si o'z only). Afterwards, we will discuss the 
magnetic field term "Hm and its effect. 

The derivation of the quantum Ising Hamiltonian 
[20; 118[ involves a slightly different interaction picture 
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we adapt the calculations from sections 3.2 and 3.4 ac- 
cordingly: The frequencies in the transformation equa- 
tion (3.14) are changed to Wm ~^ ^^s- As a result, the 
substitution e^"^™* — )- e*'"^* has to be applied to equa- 
tion (3.15) (and subsequent equations) and e^"^"* — > 1 
to equation (3.32). Hence, reads in the new interac- 
tion picture (including the expansion to first order in the 
Lamb-Dicke parameters and the RWA): 

V//(LDR)(RWA) 

N 3N 

= J2I1 i^4S^.:'°'"''«L («oi« + asaW) + h.c. 

■i^l m— 1 

(3.41) 



Figure 9. Parity measurement after the geometric phase gate 
on two radial modes of motion with two ions. The parity is de- 
fined as P :— Pii + P-f-f — (P-fi + where P^.s' denotes the 
population of the electronic state \s,s') with s, s' G {i, ti- 
lt is measured as a function of the phase ip of the analysis 
pulse i?(7r/2,7r/2 + ip) (compare figure 7). Each data point 
represents the mean of 2500 measurements. The contrast 
C — 92.2 % is determined from the fitted curve. Consider- 
ing the populations Pn + P-if > 98 % for the entangled state 
we obtain a Bell state fidelity _F > 95 %. 



However, the full Hamiltonian in the interaction picture 
now also involves 

n's = ulUsil^ = Us. (3.42) 
where the transformation is the identity, because trivially 

To gain the form of a spin-spin interaction, we apply 
a canonical transformation (compare reference 20) to the 
Hamiltonian: 



compared to section 3.2 by substituting Ho with T-Lq,: 
T~Lo — T~Lc + Hm 

= 'He + ^ fv^salnCLm - ^ hSmalnam ■ (3.39) 



=:ns 



The term Jig is added to the interaction Hamiltonian. 

To retrieve the representation of Hs in the newly de- 
fined interaction picture, 

H's — UliisU^ with := e"'^"*/'', (3.40) 



^.(LDR)(RWA) /^,(LDR)(RWA) A 



with 



N 3N 



Uc := cxp U ^ ^ 



i—l m—1 



hSr, 



and 



(3.43) 
(3.44) 



C« := i;i4%«e'^" (aoi« + aa^^) . (3.45) 

Using the calculations from section G, the transformed 
Hamiltonian reads: 



N N 3N 

^.(LDR)(RWA) +n'^ = J2J:i: T^^^ « + 



i—l j — 1 m—1 



N N 3N ^(i)^(j) U) (3) 



i—l j—1 m—1 



(3.46) 



The Hamiltonian can be expanded into a pure ® cTz 
interaction, a "bias" term with interaction, and a con- 
stant term that can be neglected. The "bias" term acts 
as a longitudinal magnetic field and leads to a deviation 
from the quantum Ising model. At first glance, this is not 
desired and it will be treated as an error in the follow- 



ing discussion. However, by including a "bias" term in a 
controlled way we could also explore an extended phase 
diagram with the longitudinal field as an additional pa- 
rameter. 

We want to stress the similarity between the spin-spin 
interaction according to equation (3.46) and the Oz ® cfz 
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interaction discussed in section 3.4. The canonical trans- 
formation has the form of a displacement operator and 
looks very similar to the first exponential function in 
equation (3.33) (except for the time dependence of the 

I 



latter). The similarity to the geometric phase term in 
equation (3.33) can be best seen comparing the time evo- 
lution operators. As the Hamiltonian in equation (3.46) 
is time-independent, the time evolution simply reads 



^.(LDR)(RWA)^^^^^) X UsitM = exp( -i^ ^ Y: 

\ i—1 j — 1 m— 1 



N N 3N r>i^}r>U)H) U) 

S S 



X 6,nit ~to)e 



5^ 



exp i ^ 5,n{t - t^)a\na„ 



(3.47) 



r 



Before we can apply the easier time evolution of equa- 
tion (3.47), in which electronic states are decoupled from 
motional states, the state vector jV') has to be trans- 
formed from the original picture to := [4 j?/;). As Uc 
depends on the electronic state, the transformation will 
in general lead to an entangled state and the canonical 
transformation can be interpreted as dressed-state pic- 
ture (electronic states "dressed" with motional states). 
As the states for the simulation of the quantum Ising 
Hamiltonian are prepared in the original (undressed) pic- 
ture, but the Hamiltonian acts in the dressed picture, an 
error is introduced into the simulation (see, for example, 
reference 119). However, as long as the effect due to 
is small (\^2^^ "qm ai j 5rn\ ^ 1), we can use the approxi- 
mation \'4')" ~ IV")- In terms of the geometric phase gate 
this corresponds to the case, when the circles in phase 
space are small and the entanglement between electronic 
and motional states can be neglected at any time. 

The same holds for the measurements of observables: 
They are performed in the original (undressed) picture, 
in which electronic states are entangled with the motional 
states, and in general a further error is introduced in the 
simulation. However, the measurement of the states is 
typically insensitive to the motional states and involves 
a projection to one of the electronic states, for example. 



As 



the ID state: p(''' 

projector does not change under the canonical transfor- 
mation and the readout of ct^ eigenstates (without any 
rotations of the bases applied beforehand) does not in- 
troduce further errors. 

The form of the "bias" term proportional to tr^ can be 
simplified in the case of a linear Paul trap with equal 
Rabi frequencies /^g*-* and equal phases cpg'' for all ions: 
The sum over j extends over the Lamb-Dicke parameters 
only. This sum is non-zero only for centre-of-mass 



J3) 
tjm 



modes, for which the t]™' additionally are independent of 
the site j. Hence, the three sums simplify to a sum over 
with constant prefactor [20]: 




(3.48) 



However, this simplification does not necessarily hold for 
two-dimensional arrays of individual traps for each ion. 

In the following, we will discuss the magnetic field 
term, which originates from a ax interaction described by 
equation (3.17) (with ai = 1 and a2 = 0). In principle, 
we have to apply the substitution e^"^"* — >■ e^'"^* due to 
the new interaction picture here, too. However, we con- 
sider a magnetic field term without motional dependence 
in the following (riim = 0, compare section 3.3) and thus 
the terms containing the motional creation/annihilation 
operators vanish: 



(3.49) 



The canonical transformation can be rewritten as 
J7c = exp ^i^^**^ (aoi^''^ +a3ai'^) j (3.50) 



with the Hermitian operator 



Wt; 



0, the and 



Sm 



(3.51) 



(3.52) 



Trivially, the commutator 



0. The canonical 



transformation of H 



/(RWA) 
M 



(see equation (3.49)) is thus 



equivalent to a transformation of the iTj!'' operator as in 
equations (B.16) and (B.17): 

^//(RWA) _ ^/(RWA)/-.t 



'■M 



M 



N 
1=1 



h.c. 



(3.53) 
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The expansion to first order in (and thus to first order 
^^^"^^ yields: 



m 77^ 



1 + 2ia3 



3JV 

E 

m— 1 



+ h.c. 



/(RWA) 
M 



(3.54) 

The magnetic field term after the canonical transforma- 



tion deviates to order fl. 



r]m as/ 0. 

This introduces an 
additional error in the simulation. If the condition 
l^g'^'Umai/Sml ^ 1 is met, it can be small or even negli- 
gible and we effectively will obtain the desired magnetic 
field term. 

To summarize, the complete Hamiltonian is obtained 
by adding equation (3.46) and equation (3.54). It consists 
of a spin-spin interaction term and a simulated magnetic 
field pointing in x-direction, which add up to the ideal 
quantum Ising Hamiltonian. Assuming a resonant inter- 
action for the simulated magnetic field (wm — i^tJ- = 0) 
and neglecting the phases (v^g ^ ~ y^^^f = 0), the quantum 
Ising part can be written 

^Qlsing •= + 'Hj 

N N N 

= ^ sWaW J^'''^^^^ ® ^^'^ ^^'^^^ 

i=l 4=1 j = l 



the pure ax interaction 'H'^^^\ 



due to from 



where 



^{i)r.{3) [i) {3) 

B<^'^:=hn^^ and := ft ^ if^il^^^I^i^a^. 

m— 1 

(3.56) 

(Note that the superscripts of Bx'^ and J*^*'^^ indicating 
the site will be omitted in the following sections, if the 
interaction strengths for all ions are equal.) In addition 
to the Ising part we obtain the following terms that lead 
to a deviation from the ideal model (constant terms are 
omitted) : 



AT N 



ns + n'E- (3.57) 



The first "bias" term is further discussed in the context 
of the experimental realization, see reference 41 and sec- 
tion 5.1. The second term can be interpreted as an 
energy offset, which cancels by applying an appropriate 
redefinition of the energy scale. As mentioned above, the 
last term leads to only a small or even negligible error 
for \il^''r]mai/Sm\ <C 1. For a more detailed discussion 



of the errors in the simulation of quantum spin Hamilto- 
nians we refer the reader to reference 119. 



0, such that 



We want to emphasized that 
the time evolution of the total quantum Ising Hamilto- 
nian T^Qising cannot be simply described by the time evo- 
lutions of Hb and Hj separately. 



4. OPERATIONS INTERPRETED FOR 
EXPERIMENTAL QS 

To realize a QS for a quantum spin Hamiltonian, we 
have to (1) simulate the spin, provide (2) its initializa- 
tion and (3) the interaction of this "spin" with a simulated 
magnetic field, (4) realize an interaction between several 
spins (spin-spin interaction), and (5) allow for efficient 
detection of the final spin state. Additional diversity for 
QS arises by the capability of precise initialization, con- 
trol and readout of the motional states. 

The mathematical derivation and description of the 
individual operations have been described in section 3. 
In this section, we explain in a simplified pictorial way 
the related generic building blocks in terms of an adia- 
batic QS of a quantum spin Hamiltonian within a linear 
chain of ions. No specific ion species or trapping con- 
cept is required. A well-suited system to illustrate the 
generic requirements and to investigate the feasibility of 
QS in ion traps is given by the quantum Ising Hamil- 
tonian (see equation (3.55)). We want to note that the 
building blocks already suffice to implement a whole fam- 
ily of quantum spin Hamiltonians. 



4.1. Simulating the Spin 

The mutual distance between the ions/spins in linear 
RF traps is typically of the order of several microme- 
tres (see figure 3). Therefore, the direct interaction be- 
tween their electronic states remains negligible, which is 
advantageous, because the related interaction strength 
could hardly be tuned or even switched off. Therefore, 
the spin- 1/2 states are implemented like qubit states (see 
section 2.3). 



4.2. Simulating the Magnetic Field 

Implementing an artificial spin allows to shape ar- 
tificial fields to implement a precisely controllable in- 
teraction and related dynamics between the "spin" and 
the "field". To simulate an effective magnetic field, the 
two electronic states ||) and |t) are coupled via electro- 
magnetic radiation (see section 2.4, operation (a)). The 
related coherent oscillation of the state population be- 
tween the two levels can be described in terms of Rabi 
hopping. In the Bloch sphere picture, the tip of the elec- 
tronic state vector rotates during one flop continuously 
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from state \\.) to |t) and vice versa. For continuous cou- 
pling this can be interpreted as the precession of a spin 
exposed to a perpendicular magnetic field. 

The rotation matrix in equation (3.30) exactly de- 
scribes this interaction with a single spin (see also sec- 
tions 3.3 and 3.5). For example, if we start with |^) 
and apply a pulsed rotation i?(7r/2, 7r/2), we will ob- 
tain an eigenstate of ct^., which is abbreviated by |— >) := 
1 / v^(l^) + It))- In tlie Bloch picture, this corresponds to 
a 90 ° rotation of the Bloch vector around the y-axis, such 
that it will point in direction of the x-axis. Continuing 
with a second identical rotation we just flip the spin to jt) 
as if we applied R{tt, 7r/2) or a 180 ° rotation around the 
y-axis, respectively. However, we can replace the second 
operation by R{tt/2,0), which corresponds to a rotation 
around the x-axis. As the state |— >) is an eigenstate of 
ax it will not be affected. 

Stroboscopic rotations have been introduced in sec- 
tion 3.5 to implement single-qubit gates for a QC. Con- 
tinuous versions of these single-qubit operations can be 
interpreted in the context of analogue QS as simulated 
magnetic field (first term of equation (3.55)). 



4.3. Simulating the Spin— Spin Interaction 

Let us first discuss the implementation of a basic spin- 
spin interaction close to the original proposal in refer- 
ence 20: Two ions are confined in a linear RF trap and 
a standing wave provides state-dependent dipole forces. 
The ions are located at same phases ((^'^'^ = 0), such that 
ions in different spin states are pulled/pushed in opposite 
directions. 

If both ions are in the same spin state, they will 
be pulled in the same direction. Hence, their mutual 
distance and mutual Coulomb energy, respectively, re- 
mains unchanged. However, if the two spins are in dif- 
ferent states, one ion will be pulled and the other one 
pushed. Their mutual distance and as a result their mu- 
tual Coulomb energy will change. This is exactly the 
essence of a spin-spin interaction, where the energy cor- 
responding to a spin state depends on the states of its 
neighbours. To interpret interactions as ferromagnetic 
or anti-ferromagnetic it is advantageous to consider the 
mutual Coulomb energy in longer chains of spins (see 
figure 10). 

The technical realization in reference 41 avoids the 
difficulties arising from standing waves and resonantly 
enhances the interaction strengths by implementing the 
spin-spin interactions with stimulated-Raman transi- 
tions as in the case of quantum gates [118] (see sec- 
tion 3.5). In a pictorial interpretation, the standing 
waves are replaced by "walking" waves and instead of 
static displacements we obtain driven oscillations of the 
ions. However, the mathematical description yields ex- 
actly the same spin model in an appropriately chosen 
frame (see section 3.6). The sign of J can additionally 
be changed by choosing a different sign for the detunings 
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Figure 10. Ion chains superimposed by standing waves pro- 
viding state-dependent forces in (a) axial direction and (b) 
radial direction. All ions are placed at the same phases of the 
standing waves, (a) If all spins are in the same state, the ions 
will all be shifted in the same direction without changing their 
mutual Coulomb energy. However, if every second spin is in 
the opposite spin state, distances between neighbouring ions 
will be alternately increased and decreased and the mutual 
Coulomb energy is increased due to its 1/d dependence. Here 
d denotes the distance between neighbouring spins. The fer- 
romagnetic order is energetically preferred, such that J < 
for this interaction, (b) A chain of ions with same spin states 
will again only be displaced and the mutual Coulomb energy 
will not change. For alternating spin states the distances be- 
tween neighbouring ions will increase, such that the mutual 
Coulomb energy will be decreased. (Note that this should 
not to be confused with the structural zigzag phase-transition 
(see figure 3), where the the radial displacements are typi- 
cally orders of magnitude larger and spin-independent.) The 
anti-ferromagnetic order is energetically preferred, such that 
J > 0. 



6m from the modes (see equation (3.56)). 



4.4. Geometric Phase Gate versus Adiabatic QS 

It might be helpful to emphasize the differences and 
similarities of the interactions in QC and analogue QS: 
To realize a phase gate operation on the radial modes of 
two qubits, as described in section 3.5, typically one or 
two isolated modes of motion are selected. The small de- 
tuning from the modes is chosen to obtain comparatively 
large interaction strengths and thus motional excitations. 
For ions being initialized in the motional ground state 
the displacements in the respective phase space(s) lead 
to an average phonon number n 1 and to a significant 
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entanglement between electronic and motional states at 
intermediate times. However, this entanglement can be 
withdrawn until the end of the gate and maximized the 
final entanglement between the qubit states only (see sec- 
tion 3.5). 

In contrast, we consider an adiabatic evolution accord- 
ing to the quantum Ising Hamiltonian in the case of 
analogue QS (see section 5.1 for the experimental pro- 
tocol). We have to make sure that the entanglement 
between electronic and motional states remains small at 
any time during the simulation (see discussion of errors 
in section 3.6). Additionally, running the simulation on 
many spins simultaneously will result in contributions 
from many motional modes simultaneously. As a result, 
a large detuning from all modes has to be chosen, such 
that the difference of the radial frequencies can be ne- 
glected and a net effect from all modes remains. Choos- 
ing the right parameters allows to simulate spin-spin in- 
teractions of different strength, different signs and even 
range of interaction [20] . 

Furthermore, a scan of the duration of the displace- 
ment pulses Td in geometric phase gates leads to a pe- 
riodic evolution from to a-^^d vice versa (see fig- 
ure 8). 

In contrast, the distinct contributions (Hb and "Hj) of 
the quantum Ising Hamiltonian are not stroboscopically 
alternated but applied simultaneously. As mentioned in 
section 3.6, the time evolution according to the quan- 
tum Ising Hamiltonian is not simply the time evolution 
according to Hb and "Hj separately. As a consequence, 
applying the spin-spin interaction for a longer duration 
and/or increasing its strength does not alter the state 
anymore. 



4.5. Note on Simulating (Virtual) Particles for QS 

Up to now we summarized how the tools developed for 
QC can be adapted, used and interpreted as tools for ana- 
logue QS. However, the toolbox for QS is substantially 
larger (see also section 5.2). (1) Phonons do not have to 
be restricted to mediate interactions in QC and QS: They 
where also proposed to simulate bosons, for example 
atoms in the Bose-Hubbard model [21] or charged parti- 
cles [120]. (2) Topological defects in the zigzag structure 
of two-dimensional Coulomb crystals (see figure 12) are 
proposed to simulate solitons ]121]. 



5. TOWARDS SIMULATING MANY-BODY 
PHYSICS 

In the first part of this section we want to assemble 
the building blocks described above to illustrate, how an 
analogue QS of a quantum spin Hamiltonian can be im- 
plemented. For this purpose, we will describe the re- 
alization of first proof-of-principle experiments on the 
quantum Ising Hamiltonian (see equation (3.55)). In the 



second part we aim to summarize, to the best of our 
knowledge, the existing proposals addressing many-body 
physics with the described and available toolbox. 



5.1. 



Proof-of-Principle Experiments on Quantum 
Spin Hamiltonians 
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Figure 11. Probability to find two spins in either of the states 
I44-) or Itt) after an adiabatic QS of the quantum Ising Hamil- 
tonian in dependence of \ J/Bx\, starting within paramagnetic 
order. The experimental protocol (top) consists of the in- 
teractions applied simultaneously including an adiabatic in- 
crease of |J| to transfer the system from the former ground 
state |— >■—>■) to the new one (bottom). We achieve a maxi- 
mal probability of P^f = Pj_j_ — (49 ± 1) % to observe one 
of the states |4^) and |tt) corresponding to a ferromagnetic 
order and define the quantum magnetization to be equal to 
Ptt + -^1-4- = (98 ± 2) %. We derive the fidelity for the entan- 
gled state l/\/2(|||) + Itt)) to approximately F = 0.88 by a 
parity measurement (compare reference 41 and section 3.5). 

First, we will describe the basic implementation of the 
experimental protocol on the axial modes for the case of 
two spins ]41], as illustrated in figure 11. Subsequently, 
we will emphasize the differences and additional informa- 
tion explored in references 43 and 42. For the details on 
the individual experimental parameters we refer to these 
references. 

For the case of two spins, the protocol has been realized 
following five steps: (1) The two ions are initialized by 
Doppler cooling, sideband cooling, and optical pumping 
(see section 2.5) in the state HI) |nsTR = 0). (2) Both 
spins are prepared by a common i?(7r/2, 7r/2) rotation 

in the o-i'' eigenstate |— I^istr = 0). (3) An effective 
magnetic field of amplitude is applied equivalent to a 
continuous i?(2i?Mi>0) rotation (see equations (3.55) and 
(3.56)). At this step, the state |— >■—>■) I^str = 0) repre- 
sents the ground state of the first term of the quantum 
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Ising Hamiltonian in equation (3.55) that can be "easily" 
prepared. Note that the rotation is shghtly ofF-resonant 
to mimic an additional interaction counteracting the 
"bias" field (see also reference 41). (4) The effective spin- 
spin interaction J is ramped up adiabatically with re- 
spect to the timescale 1/X2m defined by the simulated 
magnetic field, until | J| ^ B^. The system adiabatically 
evolves into its new ground state, which is an equal su- 
perposition of the two energetically preferred states of 
the ferromagnetic order: l/v^(||4.) -I- |tt))- (5) Finally 
both interactions are switched off. The readout of the fi- 
nal spin state is performed by state dependent detection. 
This projects the spin state to one out of the four eigen- 
states of the measurement basis (|||) , \\X) , |tl) , |tt))- 
Steps (1) to (5) are repeated many times to obtain the 
populations related to these states. 

To investigate the degree of entanglement of the fi- 
nal spin state, an additional parity measurement is per- 
formed as in the case of the geometric phase gates (see 
section 3.5). The populations of jH) and \\X) in depen- 
dence of \ J\lBx and the entanglement fidelity are sum- 
marized in figure 11. 

The experimentally observed entanglement of the fi- 
nal states confirms that the transition from paramagnetic 
to ferromagnetic order is not caused by thermal fiuctua- 
tions that drive thermal phase transitions, but by the so- 
called quantum fluctuations [39, 122] driving QPTs in the 
thermodynamic limit at zero temperature. In this pic- 
ture tunnelling processes induced by coherently cou- 
ple the degenerate states |t) and jj,) with an amplitude 
oc Bx,l\J\. For N spins the amplitude for the tunnelling 
process between = |tt • ■ • t) and \^nCi = ||| . . . |) 

is proportional to {B^I\J\)^ ^ since all TV spins must be 
flipped. In the thermodynamic limit (TV — )■ oo) the sys- 
tem is predicted to undergo a QPT ak, \ J\ = B^. At 
values I J| > Bj^ the tunnelling between l^'oot) and \i'oo\) 
is completely suppressed. In our case of a finite system 
|tf'2t) and 15^2^) remain coupled and the sharp QPT is 
smoothed into a gradual change from paramagnetic to 
ferromagnetic order (see figure 11). 

It has to be noted that the performance of such a sim- 
ulation on a large number of spins in a one-dimensional 
chain requires several technical improvements. Recently, 
the group at the University of Maryland pioneered a sub- 
stantial step for scaling by investigating the emergence of 
magnetism in the quantum Ising model using up to nine 
ions [43]. To achieve these results they mediated the in- 
teractions via the radial modes of motion [20, 116[ (see 
also section 3.5). Furthermore, they implemented the 
effective spin-spin interactions in a rotated frame using 
M0lmer-S0rensen interactions [79, 81] on robust hyper- 
fine clock states. Thereby, they do not depend on the 
phases t/j^'^ of the laser beams at the sites of the ions. To 
perform an adiabatic transition, the simulated magnetic 
field has been adiabatically turned off, while the effective 
spin-spin interaction remained constant. 

Their results allow already much more than simply in- 
creasing the number of spins: they enter a new regime of 



intriguing questions. The crossover of the quantum mag- 
netization ]41] from paramagnetic to ferromagnetic order 
is sharpening as the number of ions is increased from two 
to nine, "prefacing the expected quantum phase transi- 
tion in the thermodynamic limit" [43]. Even though the 
results can still be calculated on a classical computer, 
they provide a possibility to critically benchmark QS 
aiming for only slightly larger, but intractable systems. 

Already increasing the number of ions to three and 
adapting the individual spin-spin interactions includ- 
ing their signs allows to address spin frustration in the 
smallest possible magnetic network ]42]. Spin frustra- 
tion of the ground state can be pictorially understood in 
a two-dimensional triangular spin lattice featuring anti- 
ferromagnetic spin-spin interactions. Here, it becomes 
impossible for neighbouring ions to have pairwise op- 
posite states. Classically, two ions will adopt different 
states, while the state of the third one is undetermined. 
During an adiabatic evolution of the quantum mechan- 
ical system (starting from the paramagnetic order) na- 
ture will choose a superposition of all degenerate states, 
leading to massive entanglement in a spin-frustrated sys- 
tem. In the realization of the experiment, the three ions 
are still trapped in a one-dimensional chain. However, 
almost complete control over the amplitudes and signs 
of 7*^^'^^ is gained by coupling to particu- 

lar collective modes of motion and choosing appropriate 
detunings ]42]. 

It has to be mentioned that for an increased number of 
spins the energetic gap between ground and excited states 
further shrinks and the requirement on adiabaticity en- 
forces longer simulation durations related to a longer ex- 
posure to decohering disturbances. Still, as mentioned 
in section 1, the influence of decoherence might destroy 
the entanglement within the system, but this might not 
be relevant for the observable of interest. Here it will be 
crucial to investigate the role of the decoherence effects 
with respect to the specific analogue QS. 

With respect to digital (stroboscopic) QS it should be 
emphasized, that no quantum error correction is required 
for proof-of-principle experiments on a few ions. Promis- 
ing results of a "trotterized" version of the simulation of 
the quantum Ising Hamiltonian with two spins have been 
performed recently [123]. 



5.2. Systems Featuring Many-Body Physics 
Proposed for Analogue QS 

Condensed matter met atomic, molecular, and optical 
physics not so long ago, when trapping techniques for ul- 
tracold neutral atoms and ions allowed experimentalists 
to generate lattices and crystals, where models from solid 
state physics may be implemented. Combining the fields 
has led to a very rich interdisciplinary research activity, 
as well as to several misunderstandings between scien- 
tists looking at the same system from different points of 
view. In the particular case of trapped ion experiments. 
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the outlook for quantum simulation of many-body mod- 
els is very exciting, but some knowledge on the details of 
this physical system is required to understand both the 
limitations and the amazing possibilities of this setup. 

In the following we review many-body models that 
have the potential to be simulated with trapped ions. 
There have been several contributions both from theory 
and experiments to this research line. Most of them share 
the common feature that they are inspired by known 
models from condensed matter physics, but their imple- 
mentation with trapped ions turns out to lead to a rich 
variety of new physical phenomena, which may even re- 
quire new theoretical paradigms that go beyond the con- 
ventional ones in the solid state. The three main reasons 
for that are: (1) Trapped ion experiments are naturally 
performed in a non-equilibrium regime, whereas solid 
state physics typically deals with thermal equilibrium, 
(2) trapped ion systems may in principle be controlled 
and measured at the single particle level, (3) ion crystals 
are typically mesoscopic systems, in the sense that they 
may reach a number of particles (spins, phonons, . . . ) 
large enough to show emergent many-body physics, but 
still finite size effects are important. All those peculiari- 
ties have to be kept in mind, since they provide us with 
unique features for analogue QS. 



5.2.1. Quantum Spin Models 

Following the experimental advances in QIP, the most 
natural degree of freedom to be used for QS seems to be 
the electronic states for spins and the phonons to mediate 
their mutual interactions. However, one has to identify 
conditions where interesting phenomena arise, such as, 
for example, quantum critical phases. This has already 
lead to the promising proof-of-principle experiments dis- 
cussed above. 

A unique feature that we can exploit with trapped ions 
is the fact that the effective spin-spin interactions can be 
implemented showing a dipolar decay, J^''-') cx l/|i — jj'^. 
In the case of the Ising interaction, the cubic depen- 
dence does not change the critical universality class of the 
model, as shown, for example, by the numerical calcula- 
tions in reference 119. However, even in this case, long- 
range entanglement is induced by the long-range interac- 
tion, which is absent in conventional nearest-neighbour 
quantum Ising chains. On the other hand, when con- 
sidering other interacting schemes, like the xyz Hamilto- 
nian, the dipolar interaction may lead to the formation 
of quasi-crystalline phases of spin excitations [124]. 

Several pieces have been added to the toolbox of quan- 
tum simulation, which definitely allow us to explore 
physics beyond conventional solid state paradigms. For 
example, a theoretical proposal has been presented to 
implement models, whose ground states show topolog- 
ical features [125]. Also, methods to implement three- 
body spin-spin interactions have been designed, see ref- 
erence 126. Finally, dissipation in trapped ion systems 



has been proved to be useful to engineer quantum phases 
that arise as steady-state of dissipative processes [127]. 
The many-body physics of dissipative systems is a much 
more unexplored area than equilibrium properties, even 
for theorists. For that reason, adding dissipation to quan- 
tum magnetism opens an exciting perspective for trapped 
ions. 



5.2.2. Interacting Boson Models 

A variety of exciting quantum many-body systems may 
also be simulated by using the collective motional de- 
grees of freedom (phonons) to realize models of interact- 
ing bosons. In particular, whenever the motional cou- 
pling between ions is small compared to the trapping fre- 
quency, the phonon number is conserved and becomes 
a good quantum number to characterize the quantum 
state of the system. This principle was introduced and 
exploited in reference 21 to show that the physics of ra- 
dial modes in Coulomb chains is effectively described by 
a Bose-Hubbard model. Vibrational couplings between 
two ions, say 1 and 2, induced by the Coulomb interac- 
tion, have a typical form cx x^^''x^'^\ where x*^'-* is the 
ion displacement operator. Under the approximation of 
phonon number conservation, those terms become tun- 
nelling couplings of the form (a|a2 -l-h.c.). The same idea 
applies to quartic anharmonicities of the trap, which yield 
Hubbard interactions, (aJ^Om)^- Anharmonicities may 
be induced and controlled with optical forces, as shown 
in reference 21. This analogy between phonons and in- 
teracting bosons opens an exciting avenue of research, 
where experiments might be relevant even with a single 
ion, realizing a single anharmonic quantum oscillator. 

The ground state of those phonon-Hubbard models in 
Coulomb chains was extensively studied in reference 128, 
where it was shown that phonon Luttinger liquid phases 
may arise. Very recent experiments show indeed the tun- 
nelling of phonons between ions trapped by different po- 
tentials, realizing thus an important step towards the use 
of phonons for quantum simulation [129, 130]. Exploit- 
ing phonon tight-binding models has been also shown to 
allow us to implement models with disorder showing An- 
derson localization [131], as well as synthetic gauge po- 
tentials by using periodic driving of the trap frequencies, 
see reference 120. Dipole forces acting on ions confined 
in a microtrap array (see section 6), motional couplings 
can be controlled such that phonons experience synthetic 
gauge fields. This idea would lead to the simulation of 
magnetism in quantum lattices with trapped ions. 

5.2.3. Spin-Boson Models 

The natural convergence of the proposals presented 
above leads to the quantum simulation of spin-boson 
models. This is a paradigmatic model for quantum im- 
purities in solids, which typically describes a single spin 
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coupled to a continuous bath of harmonic oscillators with 
a power-law spectral density. Surprisingly, the coupling 
of the electronic levels of a single ion to the axial phonons 
of a Coulomb chain yields a spin-boson model with a 
quasi-ohmic spectral density [38] . The physics to be sim- 
ulated here is equivalent to some celebrated models in 
condensed matter physics, such as the Kondo effect. The 
finite size effects that are intrinsic of trapped ion sys- 
tems, turn out to yield features beyond the conventional 
physics of these models, in particular quantum revivals 
associated with the reflection of vibrational waves along 
the chain. Quite recently, it has been proposed to study 
a situation in which spins and phonons are coupled, in 
such a way that a Jaynnes-Cummings-Hubbard model 
is simulated [132]. In this model phonons follow a tight- 
binding Hamiltonian and, in addition, they are locally 
coupled to spins. The system has been shown to undergo 
a superfluid-Mott insulator QPT. 



5.2.4- Inhomogeneous Many-Body Models: Impurities and 
Topological Defects 
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Figure 12. Topological defects in two-dimensional Coulomb 
crystals (compare figure 3c for a comparable crystal with- 
out defects). Changing the experimental parameters non- 
adiabatically during a structural phase transition from a lin- 
ear chain of ions to a zigzag structure, the order within the 
crystal breaks up in domains, framed by topologically pro- 
tected defects that are suited to simulate solitons. (a) Numer- 
ical simulations for 33 ions predicting a localized topological 
defect at the position of the marked (blue) ions, (courtesy of 
Benni Reznik and Haggai Landa) (b) CCD image of 45 laser 
cooled Mg"*" ions providing clear evidence of the topological 
defect indicated by the zigzag-zagzig transition. The crystal 
contains a non- fluorescing molecular ion (MgH^) at the red 
mark, (courtesy of Ciinther Leschhorn and Steffen Kahra of 
the group at MPQ) 

The tools for QS in ion traps are not restricted to elec- 
tronic and motional degrees of freedom only. It has been 
proposed to exploit impurities in the Coulomb crystal. 
On the one hand, for example, by embedding ion(s) of 
a different species (different mass) into the crystal and 
taking advantage of the altered spectrum of the modes 
and scattering of phonons [133] and the option to include 
larger simulated spins {S > 1/2) ]134]. On the other 
hand, by creating localized topological defects within the 
more dimensional structure of the Coulomb crystal (see 
figure 12). In reference 121 it was suggested to induce 
a structural phase transition from a linear chain of ions 
(see figure 3b) to a zigzag structure (see figure 3c), for 



example, by lowering the radial confinement. Changing 
the parameters in a non-adiabatic way (fast compared to 
the phonons mediating information within the crystal) 
should cause independent domains of "zigzag" and "za- 
gzig" structure, respectively. At their clash, topological 
defects were predicted and have recently been observed 
(see figure 12). The number of the created defects should 
scale according to the Kibble-Zurek prediction ]135-137]. 
The defects themselves can be interpreted as solitons 
]121]. Solitons are defined as localized solutions of nonlin- 
ear systems, which depend essentially on the nonlinearity. 
Such solitons have a unique spectrum of frequencies with 
modes which are localized to the soliton and whose fre- 
quency is separated by a gap from the other phonons. A 
quantum mechanical time evolution of these modes was 
calculated numerically and it is expected to remain coher- 
ent for hundreds of oscillations ]138]. QS could allow to 
explore their potential applications for QIP [121[ as well 
as the quantum behaviour of these "objects" themselves. 
Solitons appear in all branches of the natural sciences and 
have been extensively investigated in solid state systems 
[139[. Among others, classical solitons were observed in 
waveguide arrays [140, 141[ and Bose-Einstein conden- 
sates [142[, where they are mean field solutions. Discrete 
solitons were investigated in the Frenkel-Kontorova (FK) 
model [143, 144[, which describes chains of coupled par- 
ticles interacting with a local nonlinear potential. In a 
different realization, a variant of the FK model can also 
be realized in the ion trap by adding an optical lattice to 
a linear chain [145-147[. 

The important requirement to address any of these in- 
triguing models will be to increase the number of ions 
and the dimensionality of the system. Trapping ions in 
two-dimensional arrays would allow to study hard-core 
boson phases, showing the effect of frustration, quantum 
spin liquid phases, and quantum states with chiral order- 
ing [148[. Two approaches for scaling will be described 
in more detail in the following sections. 



6. SCALING ANALOGUE QS IN ARRAYS OF 
RF SURFACE ELECTRODE TRAPS 

One possible way to overcome the limitations on scal- 
ability of trapped ions in a common potential well (see 
section 2.2) is to store them in an array of individual RF 
traps. 



6.1. One-Dimensional RF Surface Electrode Traps 

Conventional RF traps with their three-dimensional 
geometry of electrodes (see figure 2) individually fabri- 
cated with conventional machining were unique "master- 
pieces" with unique characteristics. 

In 2005 and 2006, the group at NIST pioneered the 
miniaturization of RF traps by projecting the electrodes 
onto a surface [67, 68[ (see figure 13a.), very similar to 
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Figure 13. Schematic to illustrate the projection of the elec- 
trodes of the RF (yellow) and DC (shaded) electrodes on a 
surface as a way to scale towards two-dimensional arrays of 
ions. The black crosses indicate the positions of the minima of 
the pseudopotentials. (a) Cross section of the electrodes of a 
conventional linear RF trap with three-dimensional geometry 
and the electrode structures projected onto a surface. The 
dashed arrows point at the new location of the electrodes, the 
white areas represent isolating gaps, (b) Cross section (upper 
part) and top view (lower part) of the stripe electrodes. It 
has been proposed to concatenate several of linear RF surface 
electrode traps depicted in (a) as a basic unit to span a two- 
dimensional array of ion [155] (red and blue disks representing 
ions in opposite spin states). For sufficiently small mutual ion 
distances and decoherence rates of the ions, this approach is 
intended to scale analogue QS. 



chip traps for neutral atoms [149]. Introducing photo- 
lithographic techniques for the trap fabrication opened 
up exceptional precision and the production of small se- 
ries of identical traps, see for example references 150 and 
151. Within these linear RF surface electrode traps, mo- 
tional ground state cooling was achieved at a height of 
the ion over the electrode surface of ft- = 40 pm and with 
a comparatively small motional heating rate of the order 
of 1 quantum/ms [68[. 

Motional heating rates scale with ^ ft""* [152[. The 
exact heating mechanisms are not yet fully understood 
and the groups at NIST, in Berkeley, at MIT and oth- 
ers currently put a lot of effort in further investigations. 
However, the groups at MIT [153[, NIST, and the Uni- 
versity of Maryland demonstrated a significant reduction 
of the heating rates in cryogenic (surface electrode) traps 
for QC purposes (see also reference 44). The inverse of 
these heating rates is long compared to typical opera- 
tional durations of a QC of tens of microseconds (see 
also section 3.5). 

For scaling towards a universal QC it might suffice 
to interconnect linear ion traps via junctions on a two- 
dimensional surface to a network of one-dimensional 
traps [154[, realizing the "multiplex ion trap architecture" 
[69[. That is, ions are proposed to be shuttled between 
processor and memory traps only interacting in the pro- 
cessor traps. This would allow to subdivide the large 
total number of ions into small groups in many indi- 
vidual traps and to reduce the local requirements to a 
technically manageable effort. One-dimensional RF sur- 
face electrode traps with more than 150 individual DC 
electrodes and several junctions have been realized [154[, 



allowing to shuttle ions at moderate heating rates. 

In addition the opportunity arose to deliver identical 
traps to different groups. One example is the linear RF 
surface electrode trap (denoted by "Sandia Linear Trap" 
in the following) [150[, which was designed by the groups 
of Oxford, Innsbruck and Sandia National Laboratories. 
The latter fabricated a small series of identical replicas. 
The traps have been tested in several laboratories and the 
individually measured trapping parameters are in good 
agreement with the design values. Publications are in 
preparation by the groups at Oxford and Sandia (see also 
reference 150). 

It has to be emphasized that pursuing the multiplex 
approach for scaling universal QC is not applicable to 
the proposed analogue QS, where the ensemble of spins 
is supposed to evolve uniformly as a whole. 



6.2. Optimized Two-Dimensional Arrays of RF 
Surface Electrode Traps 

Shortly after the invention of RF surface electrode 
traps it was proposed to concatenate linear traps suf- 
ficiently close, such that the ions experience mutual 
Coulomb interaction in two dimensions [155] (see fig- 
ure 13b). However, for a real two-dimensional lattice at 
sufficiently small and uniform ion distances of d < 40 pm 
in two dimensions, this proposal requires the ions to ap- 
proach the disturbing surface to h < d/2 ^ 20]im [67[. 

R. Schmied et al. implemented a method to calculate 
the global optimum of the electrode shapes for arbitrary 
trap locations and curvatures (originally only for periodic 
boundary conditions) [156]. The gaps between neigh- 
bouring electrodes were neglected. The authors exem- 
plary optimized a trap array with comparatively stiff hor- 
izontal confinement. 

The idea of optimizing electrode structures can also be 
used for designing traps for analogue QSs with partially 
converse requirements. In a collaboration of R. Schmied, 
NIST, Sandia National Laboratories, and us, such a sur- 
face electrode trap has been designed and is currently in 
fabrication. The trap will provide three trapping zones 
arranged in a triangle (similar to figure 14) and is in- 
tended as a first step towards larger arrays of ions. For 
this purpose, the optimization method was extended to 
finite-sized traps. 

It has to be emphasized that there are currently several 
proposals and approaches for arrays of surface electrode 
traps mainly for QC. Groups in Berkeley and Innsbruck 
aim at trap arrays with individually controlled RF elec- 
trodes. They have the advantage of selectively lower- 
able trap frequencies for individual traps and thus in- 
creasable interaction strengths between ions in different 
traps, while especially the height of the ions above the 
surface can be larger and the trap depth of other traps 
can be sustained [157[ (see also discussion in the following 
subsections). This approach can in principal be extended 
to quasi micromotion-free shuttling of ions in arrays of 
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RF traps [158] at the expense of a precise control of the 
RF voltage for each RF electrode. Another proposal sug- 
gests individual coils to be included for each trap to allow 
for laser-less interactions mainly for QS [159]. Different 
trap geometries specifically for QSs are designed in the 
group in Sussex [160]. Arrays of Penning traps with sur- 
face electrodes are advanced by the groups at Imperial 
College [161] and the University of Mainz [162[. 

In the following subsections we discuss the optimiza- 
tion goals for a surface electrode trap for an analogue QS, 
their implications and the perspectives for scalability of 
this approach. 



6.2.1. Maximization of Interaction Strengths 

The crucial prerequisite for QSs is to maximize the 
interaction strength (see section 3), however, opposed to 
QC between individually trapped ions, while still outrun- 
ning decoherence rates. The increased mutual ion dis- 
tances in arrays of individual traps substantially reduces 
the strength of the effective spin-spin interaction. It has 
be taken into account that the conditional forces have a 
limited strength, for example, because the laser power is 
limited or the assumptions in the theoretical model im- 
pose constraints as for the quantum Ising Hamiltonian 
(see section 3.6). However, a reduced stiffness of the in- 
dividual potentials compared to the example [156[ (trap 
frequencies on the order of 27r x 20 MHz) results in larger 
displacements of the ions by the same forces. This is re- 
lated to an increased mutual Coulomb energy and thus 
larger interaction strengths. Still, a lower bound for the 
trap frequencies (on the order of 27r x 1 MHz for Mg^) 
is imposed by the constraints for efficient ground state 
cooling. 



6.2.2. Mimrnization of Decoherence 

The ions will inevitably approach the disturbing elec- 
trode surfaces, if the distance between the individual 
traps is reduced. We now reinvest the reduced require- 
ments on the stiffness of the horizontal confinement to 
increase the height of the ions above the surface h keep- 
ing the mutual ion distances d constant. Some results 
for the scenario of a basic triangular lattice are depicted 
in figure 14, which demonstrates the adapted shape of 
the electrodes due to different optimization goals. Note 
that the influence from electrodes of neighbouring traps 
increases for an increased height h. The optimization al- 
lows for an increase of the height by more then a factor of 
two. still maintaining realistic trapping parameters (see 
below). Hence, the related motional heating rates (in 
units of energy per time) are expected to be reduced by 
more than an order of magnitude. In addition, the in- 
creased h should help to protect the electrodes from the 
high intensity of the laser beams parallel to the electrode 
surfaces. 
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Figure 14. Iliustration of the optimization results for the elec- 
trode structure for a basic triangular lattice in dependence of 
the height of the ions above the traps at constant mutual 
ion distance. The white gap isolates between RF and DC 
patches. The three red disks symbolize three ions at a con- 
stant distance of d = 40 pm, hovering above the surface in a 
height of (a) h = 30 pm, (b) h = 40 pm, and (c) h — 50 pm. 
(courtesy of Roman Schmied) 

We additionally include required isolating gaps be- 
tween electrodes into subsequent simulations to deduce 
deviations in the resulting trapping potential [163] (see 
figure 14). The influence of the gaps turned out to be 
negligible for the example shown in flgure 14, however, 
for further miniaturized traps these influences will gain 
of importance due to technically limited gap sizes. 

6.2.3. Maximization of the Lifetime of Trapped Ions 

The reduced frequencies and increased height above 
the surface come at the price of a reduced trap depth. 
First, sufficiently deep potentials have to be provided 
to assure adequate loading rates out of thermal atomic 
beams, preferably via efficient photo-ionization [164, 
165[. Second, sufficient lifetimes for many ions within the 
potentials of scaled traps have to be achieved. Currently, 
the average lifetime in a room temperature surface elec- 
trode trap exceeds one hour (for the Sandia Linear Trap 
operated in our laboratory). 

Deeper trapping potentials for surface traps were al- 
ready achieved by a conductive mesh with controlled volt- 
age (85 % transmittance) few millimetres above the elec- 
trode surface [166]. It has also been successfully tested 
for the Sandia Linear Trap. The mesh shields the ions 
from charges on the window and provides a wavelength- 
independent alternative to a conductive coating (see, for 
example reference 167). 

6.2.4. Control of the Symmetry of Interaction 

We additionally gain control over the individual orien- 
tations of trap axes or the relative orientations of axes of 
different traps, respectively (see flgure 15). This allows 
to shape the interaction for a given direction of motional 
excitation between ions in different traps (see section 3). 
It also allows for cooling of all spatial degrees of freedom 
with laser beams, which have to propagate parallel to the 
trap surface to minimize scattering off the surface. We 
can rotate the individual trap axes from pointing towards 



25 




Figure 15. Electrode structures of basic triangular lattices 
with different orientation and tilt of the principal axes. Red 
disks symbolize ions trapped in the potential minima for pa- 
rameters comparable to those in figure 14b. (a) One principal 
axis points in vertical direction with the X and Y axis lying in 
the horizontal plane of the electrodes. The X-principal axis 
of each trap points towards the centre of the triangle, (b) 
The respective principal axes of all traps point in the same 
direction and additionally the Z-axis is tilted with respect 
to the surface by more than 10 ° , which results in a different 
symmetry of the electrodes. The tilt of the Z-ajds is essen- 
tial to reach all spatial degrees of freedom with laser beams 
restricted to a plane parallel with respect to the electrodes, 
(courtesy of Roman Schmied) 



the centre of the structure (see figure 15a) into a paral- 
lel alignment and additionally include the required tilt 
of the vertical {Z) axes, which will result in a different 
symmetry of the electrodes (see figure 15b). 



ions of 2tt x 2 MHz, the interaction strengths exceed by 
far 27rfi. x 1 kHz. 

In a diflterent approach, we could think of using the 
motional degrees of freedom for QSs. This scheme would 
have the advantage that bare motional couplings are 
already in the 5 kHz regime. In that sense, they are 
stronger than effective spin-spin interactions, since the 
latter are slowed down with respect to the original mo- 
tional couplings by the requirement of adiabaticity. A 
recent theoretical proposal by some of us has shown that 
by using periodic modulations of the trapping frequen- 
cies, some phenomena from solid state physics may be 
simulated, such as photon assisted tunnelling [120] (see 
section 5.2). 




Figure 16. Electrode structures for RF surface electrode traps 
scaled for analogue QSs. Black dots symbolize the RF min- 
ima, red lines serve as guide to the eye to emphasize the lattice 
structures for (a) three, (b) twelve, and (c) an infinite number 
of ions/spins, the latter both considering the parallel orienta- 
tion and the tilt of the principal axis (see figure 15). (courtesy 
of Roman Schmied) 



6.2.5. Control of the Potential in Individual Traps 



First, splitting DC electrodes into several separately 
controllable segments allows for the individual compen- 
sation of displacements of the ions from the minima of the 
pseudopotential due to stray fields and space charge ef- 
fects (compare section 2.2 and see references 168 and 167 
for schemes of micromotion compensation). Second, for 
further scaling, these electrodes can be used to compen- 
sate boundary effects. Due to larger number of inner ions, 
outer ones would be shifted to larger mutual distances. 
The further increased density of electrodes on the surface 
requires their connections in a multilayer structure with 
vertical wiring (vias) [150, 151[. 



6.2.6. Estimation of Parameters 

We estimate the strength of simulated spin-spin in- 
teractions for the case of Mg^ ions in such devices with 
currently available laser equipment. We assume an typi- 
cal laser power of 400 mW (max. 600 mW are available) 
at 280 nm from an all solid state laser source [169[. We 
further assume the beam to have a cylindrical profile with 
waists of 10 Jim x 100 jim and an electrode structure as 
depicted in figure 15b {d = h ~ 40 pm). For a trap depth 
of 100 meV and a minimal oscillation frequency of the 



6.3. Perspectives of Our Approach 

As depicted in figure 15b, in a first step three ions will 
reside on the vertices of a triangle and the interaction be- 
tween the spins can be simulated as in references 41 and 
42 (compare section 3.6) or reference 120 (compare sec- 
tion 5.2). The above parameter estimates should already 
suffice for proof-of-principle experiments and mesoscopi- 
cally scaled QS. Motional modes in two-dimensional trap 
arrays will behave similarly to radial modes in linear RF 
traps for all three dimensions [20, 118] and the effective 
spin-spin interaction will prefer anti-ferromagnetic order 
for far, red detuning from all modes. Thus, the systems 
should give us the possibility to study spin frustrations in 
a spatial, triangular configuration (see also section 5.1). 

Based on the results of these investigations further 
scaling of the surface trap architecture to large-scale (tri- 
angular) lattices of tens or even hundreds of spins might 
be pursued (see figure 16). Besides of the optimization 
of the trapping parameters, further technical difficulties 
have to be considered: 

Decoherence due to motional heating as a result of 
the vicinity to the electrode surfaces could be mitigated 
within a cryogenic setup [44, 153[. The reduced vacuum 
pressure could additionally help to increase the lifetime of 
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Mg"*", which is currently hmited by photo-chemical reac- 
tions with hydrogen (mostly H2 -l- Mg+* — > MgH+ + H* ) 
and collisions with heavy components of the rest gas. 
The reaction can also be inverted by pulsed laser beams 
[170]. However, scaling the system to tens or hundreds 
of ions will still require frequent and efficient reloading. 
Increasing the loading efficiency and preserving the vac- 
uum conditions could be achieved by photo-ionizing cold 
atoms from a magneto-optical trap (MOT) [171]. 

Currently, the available laser power should not im- 
pose any restrictions on the realization of systems of 
few tens of ions (see reference 169 providing currently 
up to 600 mW) . Higher laser powers for magnesium are 
in reach [172-174[ and could allow for even larger arrays 
of simultaneously coupled ions. Besides of that, efforts 
in optics, for example arrays of lenses [175, 176[, fibres 
integrated into the trap [177[, or integrated mirrors [178- 
180[, could provide individual addressing and high light 
intensities at the position of the ions. To further mitigate 
the problem of scattered light from surfaces, one could 
think of realizing surface traps on partially transparent 
substrates [68]. Alternatively, laser- less coupling could 
be used as mentioned in section 2.4 ]82, 83, 85-87, 159]. 

Last but not least, it still has to be identified how 
to measure observables that permit the verification of 
frustration effects without the need for full (exponentially 
complex) state tomography. 



7. SCALING QS BASED ON IONS IN OPTICAL 
LATTICES 

Some groups aim to merge the two fields of QS based 
on ions in RF traps and atoms confined in optical lattices. 
It had already been proposed to combine Coulomb crys- 
tals in a harmonic confinement of a common RF trap of 
three-dimensional geometry with (commensurate) optical 
lattices to shape anharmonic trapping potentials provid- 
ing new possibilities to simulate interactions [148]. An- 
other proposal deals with the simulation of the Frenkel- 
Kontorova model using a standing wave aligned with the 
trap axis ]147]. C. KoUath et al. suggested to exploit a 
trapped ion to coherently couple (like a scanning micro- 
scope) to the atoms confined in an optical lattice [181]. 

Optical ion trapping was realized with a single Mg'*' 
ion trapped in a dipole trap [74]. We can now dream 
of spanning an array of ions (even simultaneously with 
neutral atoms) within an optical lattice. It has to be em- 
phasized that the smaller trap depth of optical traps (see 
figure 1) renders it highly unlikely that optically trapping 
charged atoms will allow to outperform the achievable 
trapping parameters or coherence times of both, ions in 
RF traps and of optically trapped neutral atoms. How- 
ever, in our opinion, this is not required. The advantage 
of equally and closely spaced traps might be combined 
with individual addressability and, most important, the 
long-range interaction provided by Coulomb forces be- 
tween the ions. 



(a) (b) (c) 




Figure 17. Illustration of new options for analogue QS based 
on ions (red and blue) and atoms (green) in optical potentials 
(black lines as guide to the eye), (a) Ions populate an op- 
tical lattice on well separated sites. The Coulomb force still 
provides a large strength of dipolar (long-ranging) interaction 
allowing for analogue QS on many-body effects, similar to the 
proposed approach in arrays of RF surface electrode traps 
(see section 6). (b) An ion and atoms populate a common 
optical lattice and, for example, share the charge via tun- 
nelling electrons, (c) An ion could be cooled sympathetically 
by cold atoms (for example, a BEC indicated as green ellipse) 
[183, 192]. Since the micromotion of the ion and the related 
differential motion between atoms and ion becomes negligible 
in the common optical trap [54], deep equilibrium tempera- 
tures are predicted to be achievable, down to a regime, where 
ultra-cold chemistry might dominate the collisions. 

In the following, we will first describe how trapping 
of an ion in a dipole trap was achieved. Still facing a 
huge variety of challenges, the new possibilities will be 
discussed afterwards. 



7.1. Trapping of an Ion in a Dipole Trap 

The procedure used in reference 74 to load a magne- 
sium ion (^'*Mg~'') into a dipole trap consists of the fol- 
lowing steps: An atom is photo-ionized out of a thermal 
beam and trapped and Doppler cooled in a conventional 
RF trap. Next, stray electric fields are minimized at the 
site of the ion using the ion as a sensor. Then a Gaus- 
sian laser beam providing the dipole trap is focused onto 
the ion and the RF drive of the RF trap is switched 
off. From that time on, the ion is confined in the dipole 
trap in the directions perpendicular to the beam propa- 
gation. The depth of the dipole trap potential amounts 
to Uo ~ 27rfi X 800 MHz or [/q ~ x 38 mK, respectively, 
and the detuning of the dipole trap beam from the rel- 
evant transition (S1/2 O- P3/2) to Z\ « —6600/^, where 
1/r determines the lifetime of the IP3/2) state. Static 
electric fields provide the confinement in the direction of 
beam propagation. After a few milliseconds the RF drive 
is switched on again and the presence of the ion can be 
verified via its fluorescence during Doppler cooling. 

For the given parameters a half-life of approximately 
2.5 ms is achieved. This value is in very good agreement 
with the theoretical predictions, assuming exclusively the 
heating process related to off-resonant scattering of the 
trapping light by the ion, the so-called recoil heating. It 
can be concluded that the heating and subsequent loss of 
ions from the optical potential is not dominated by heat- 
ing effects related to the charge of the ion, for example, 
due to the vicinity of electrodes or fluctuating stray elec- 
tric fields. Thus, state of the art techniques for neutral 
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atoms should allow to effectively enhance the lifetime and 
coherence times [182]. 



7.2. Lifetime and Coherence Times of Optically 
Trapped Ions 

We aim to increase the lifetime by cooling the ion in the 
dipole trap. Due to the large AC Stark shift and its large 
position dependency, simple Doppler cooling within the 
existing setup is challenging. Possibilities of cooling the 
ions directly towards the ground state of motion within 
the dipole potential are currently investigated theoreti- 
cally and experimentally. 

An alternative approach suggests to use cold atoms or 
even a BEC to sympathetically cool ions [183]. On longer 
timescales the approach of cavity assisted cooling of ions 
in conventional RF traps reported in reference 184 might 
also provide long lifetimes without affecting the electronic 
state of the ion. 

Currently, the coherence time of the electronic state of 
the ion is limited to few microseconds due to the high 
spontaneous emission rate. If longer coherence times 
are required (which is not necessarily the case for ev- 
ery scenario), they can be achieved in two ways: (1) As 
for two-photon stimulated-Raman transitions the spon- 
taneous emission rate can be reduced by increasing the 
detuning. A larger beam intensity could sustain the po- 
tential depth. (2) Another option would be to work with 
blue detuned light, where the potential depths can re- 
main identical, however, the ions seek for low intensity 
and exhibit less spontaneous emission. 



7.3. Towards Ions and Atoms in a Common Optical 
Lattice 

It has still to be demonstrated that one or several ions 
can be confined within one- or more-dimensional optical 
lattices. Regarding the currently available parameters, 
we expect a mutual ion distance of the order of a few mi- 
crometres only, which corresponds to one ion at approx- 
imately every 50th lattice site. Therefore the mutual ion 
distance could remain comparable to the distance within 
one (tight) common potential of a RF trap and smaller 
than the currently envisioned distances between neigh- 
bouring traps in the RF surface electrode trap approach 
(see section 6). 

Since the photo-ionization scheme applied so far ion- 
izes out of a thermal beam of magnesium atoms, the av- 
erage kinetic energy of the atoms is much larger than 
the depth of the optical potential and, in addition, the 
local vacuum is severely affected. The loading efficiency 
for RF traps could be largely enhanced by ionizing Mg 
atoms from a magneto-optical trap (MOT) ]185], which 
would also allow to directly load atoms into an optical 
trap. In addition, after loading neutral atoms into the 
lattice, some of them could be photo-ionized on site. 



Ions and atoms confined in a common optical lat- 
tice could offer an approach to exploit the physics of 
charge transfer reactions. This might allow for a com- 
plete new class of QS, for example, of solid-state systems, 
where atoms in a completely occupied lattice (at an ini- 
tially small density of ions) share electrons by tunnelling 
causing highly entangled states of the compound system 
with most interesting quantum dynamics governed by the 
Bose-Hubbard Hamiltonian [186, 187]. 



8. CONCLUSIONS 

In the last few years the basic building blocks for a 
scalable architecture of a quantum information proces- 
sor (QC) with trapped ions have been demonstrated for 
a few qubits. Additionally, a large variety of new tech- 
niques is already tested that might considerably extend 
the available toolbox. For example, interactions based 
on magnetic field gradients and RF fields, fibre-coupled 
optical support on chip or economically and technologi- 
cally facilitated cryogenic environments. Despite the fact 
that it will be a non-trivial challenge to scale the system 
to approximately 10^ qubits, no fundamental limitations 
can be identified so far. 

On a shorter timescale, intriguing problems might be 
studied by realizing analogue quantum simulators (QS), 
by far exceeding the capabilities of classical computers. 
They can be based on similar techniques as a potential 
QC, but with less severe constraints on the fidelity of 
operations and the number of required ions. 

Currently, available operational fidelities are predicted 
to allow for studying many-body physics, for example in 
systems described by quantum spin Hamiltonians, the 
Bose-Hubbard and the spin-boson models. First proof- 
of-principle experiments simulating Ising type interac- 
tions with few ions were already successfully demon- 
strated. 

The required increase of the the number of ions and the 
accessible dimensions is proposed within two-dimensional 
arrays of RF surface electrode traps. However, the ap- 
proach is still at the level of proof-of-principle experi- 
ments and further challenges might arise during its devel- 
opment. Alternative approaches include Penning traps or 
optical lattices. 

Even though the enthusiasm within this quickly grow- 
ing field seems to be justified, it has to be emphasized 
that efficient analogue quantum simulators still require 
more than simply scaling. Examples of other important 
challenges are: (1) To investigate carefully the influence 
of different sources of decoherence on the fidelity of the 
simulation. Thus, it must be distinguished for the dedi- 
cated application, which decoherence the simulation will 
be robust against, which decoherence can be considered 
in the simulation and which decoherence is even essential 
to be included. (2) To identify possibilities to cross-check 
the validity of the output or to benchmark it against 
other QS approaches, as soon as the achieved output is 
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not accessible with a classical computer anymore. 

In the future it might be beneficial to combine ad- 
vantages of several systems for hybrid QS. On longer 
timescales, the experiences gained by developing an ana- 
logue QS based on trapped ions might culminate in ap- 
proaches incorporating solid-state devices that might al- 
low for "easier" scaling. With the realization of a univer- 
sal QC, universal QS will also become accessible. 
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Appendix A: Normal Modes and Frequencies 

This section describes a more general derivation of 
the normal modes and frequencies compared to the one- 
dimensional treatment for the linear Paul trap as in ref- 
erence 188. The equations are extended to three dimen- 
sions and an arbitrary trap potential, as long as the po- 
tential at the equilibrium position x^"^ of each ion can be 
well approximated by a harmonic potential. 

The position of the i-th ion is expressed in the Carte- 
sian coordinates of the laboratory frame 



^0 



(A.l) 



The decomposition into the equilibrium position Xq and 



displacements x^^^ yields 
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(A.2) 



= {xo,t + Xt) ex + {xo^i+N + Xi+N) eV 

+ (2:o.j+2W + Xi+2N) ez- (A. 3) 

The Lagrangian for N ions takes the form 
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where M denotes the mass of an ion, the index of the 
partial derivatives signifies its evaluation at the equilib- 
rium positions and V denotes the potential consisting of 
the trap potential Vq and the Coulomb potentials of the 
ions: 



V = Vn 
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(A.5) 



Here, Q denotes the charge and Eq the electric constant. 

For practical purposes the trap potential can be ex- 
pressed by the harmonic terms corresponding to each ion: 



N 3 



i=i j=i 



(A.6) 

(i) 

Here, ui^ denotes the j-th frequency of the harmonically 
approximated potential of the i-th ion, (f^ ^ the unity vec- 
tor of the principle axis corresponding to and p^^^ 
the position of the local minimum of the potential for 
the i-th ion. Note that the frequencies, the vectors of the 
principle axes, and the minima of the potentials become 
equal for all ions in the special case of a linear Paul trap. 

The eigenvalues of the Hessian A := (a^;) (see equa- 
tion (A.4)) yield the squares of the frequencies ujm of the 
normal modes and its eigenvectors bm determine the ions' 
motion of the m-th mode: 

qni ^b„i- X with x := {xi,..., x^n) ■ (A. 7) 
With the abbreviations q := (gi, . . . , gsTv) and B := 
(Si , . . . , b^N^ , where the shall be understood as rows 



of _B, we can express the relation in a more compact way: 



q = Bx 4^ X ^ B q. 



(A., 



Typically, the eigenvalues and eigenvectors of A have to 
be determined numerically. 



Appendix B: Transformations of Pauli Operators 

The definitions of the Pauli operators is repeated here 
to avoid confusions concerning their normalization. 
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The Pauli operators are defined by: 



Hence, the Pauh operators in the interaction picture read 
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The Pauh operators obey the relations 
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A more convenient notation in some contexts is 
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with the normahzation as in reference 44. They fulfil the 
following relations 



(B.7) 
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The transformations of the Pauli operators into the 
interaction picture involve terms of the form: 

:= e''=*^<7,e-''^'^^ (B.9) 

The transformation leaves 1 and ct^ unchanged. The non- 
trivial cases i = X and i ~ y can be calculated using the 
Baker-Campbell-Hausdorff formula 
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with B = — IkCTz, a — (T,/„ — , 1 1 LA. 

m_ := 2i. 

The commutators are given by 
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The operators (t+ and (T_ transform: 
e''^*'<T+e-'''*^ =e'2'^CT+ 



(B.16) 
(B.17) 



Appendix C: Transformations of Motional Operators 

The creation operator a and the annihilation operator 
fulfil the relations: 

a|?i) = (C.l) 

|n) = V^m|n+ 1) (C.2) 

[a,a^=i (C.3) 

The transformation of the Hamiltonians into the inter- 
action picture requires the knowledge of the transforma- 
tion of e'«("+"^). 

It can be performed using the special case of the Baker- 
Campbell-Hausdorff formula from section B again. The 
commutators appearing in the formula are 

[a + a^ -iXa^'a] = -iX (a - (C.4) 

[a + a^, — iAa^a] = — iA [a — a^, — iAa^a] 

= (-iA)2 (a + flt) (C.5) 

[a + a\-iXa^a] ^"^ = (-iA)" (a + (-l)"a^) (C.6) 

(C.7) 

Hence, the full transformation reads 
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From this relation we can immediately derive 
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by making use of the unitarity of e''*"* °. We obtain for 
the transformation: 



{a + a' ) c 



t^"^-iAdtd 



^M^[ae-- + ate-]' 



exp (i^ [i 



]) 



Appendix D: Matrix Elements of Displacement 
Operator 

The following derivation is based on reference 89. A 
difference is that we do not restrict the displacement to 
purely imaginary A in the following. A similar derivation 
can be found in appendix B of reference 105. 

The simple Baker-Campbell-Hausdorff formula 



(D.l) 



holds for [A, [A, B]] = [B, [B, A]] = 0. 
We obtain for the annihilation operator: 
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[O else 

Using the above form of the Baker-Campbell-Hausdorff 
formula, we can rewrite the displacement operator as 



D{X) = e^"*-^*'^ = c-l^l'/^e^^^e-^*'^. 
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where lI^\x) denotes the generalized Laguerre polyno- 
mials [189]. Analogously, we obtain for n' < n: 



I n'l (n-Tl') /I \ 



(D.6) 

For values A = ir^e'"' with ?; G M, we can write equa- 
tions (D.5) and (D.6) as 
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where ri< := min(n',n) and n> := max(n',ri). 



Appendix E: System of Differential Equations of 
Rabi Problem 



The Rabi problem consists of the following system of 
differential equations: 



C2 = Ae-'^^'ci 
ci = -A*e'"*C2 

It can be solved by differentiating with respect to t 

62 = Ae-''^*ci - iwAe"''^* 
ci 

and inserting equation (E.l): 
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C2 = -ia;c2 - |A| C2 
Cl = iijJCi — |ApCi. 
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Using the ansatz q = aic"^'* we obtain the character- 
istic equations 



—K2 " '^1^2 ~ |Ap 
—Ki ~ —ujKi — |Ap, 

which have the solutions 
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Here, we introduced the abbreviation k' := \J ^ + 
The solutions of equations (E.3) and (E.4) read 
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Inserting them into the original system of differential 
equations equation (E.l), we obtain the following rela- 
tions for the constants ai^±: 



The commutator of two of these terms will trivially 
vanish for all i, j and all times t' , t", if both terms belong 
to different modes m 7^ n: 



iK2,±a2,± = Xai^± 
iKi,±ai,± = — A*a2,± 
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We replace ai_± by 02, ± using equation (E.ll) and obtain 



n';^{t'),nlP{t")\ =0 for m^n. 



(F.2) 



However, for m = n, the commutators do not vanish. 
Using the relation 



C2 = la2,+e'"'* + a2,-e-'"'*U-'"*/' 
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ci = I /i+a2,+c'"'' + fi-a2,-c-'"''') c'"*/2^ (E.14) 



where we introduced the (temporary) abbreviation fj.± := 

iK2,±/A. 

The constants 02, ± can now be expressed in terms of 
the initial values C20 '-^ C2{t = 0) and cio ci(< = 0). 
Setting f = in equations (E.13) and (E.14), we obtain 
a system of linear equations with the solutions 
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= 2isin(A - A') i (F.3) 

yields 

^ 2i^i^^li> sin (<5„,(t' - n (^W - ^t^))) . (F.4) 

The time evolution operator can be calculated using 
a Magnus expansion [190, 191]. As commutators with 
higher "nesting level" trivially vanish, the expansion sim- 
plifies to: 



Hence, 



C2 = 



Cl 



{H-C2o - cio)e''''* - {n+C2o - cio)e /2 



Ai- - Ai+ 



//+(^-C2o - cio)e"' * - ^J.-{^J,+C2o - cio)e 



(E.17) 



-\k' t 



jt/2 



and by expressing in terms of uj, k' , and A; 

C2{t) - (^cos(k'A) + |-^sin(K't)^ c-''^*/2^2(0) 

+ Asin(K'i)e-''^*/2ci(0) 
ci(i) = -^sin(«;'t)e'"*/2c2(0) 



+ ( cos(k'A) - -^sin(K'i) ) e'^'/^ciiO). 
\ 2 k' ' 



(E.18) 



(E.19) 



(E.20) 



Appendix F: Time Evolution Operator 

The calculation of the time evolution operator involves 
terms of the form 



(F.l) 



where G M and the total Hamiltonian reads 'H{t) —: 
Z^ili Z^r^i '^m (*)• (More generally, the constants 



represent Hermitian operators with 
i,j,m,n.) 



ov 



U{t,to) = e^p \ -- / dt"H{t' 



1 



to 

rt nt' 
•/to Jto 



(F.5) 



n{t'),n{t") 



The single integrals yield 



t (0 
dt'H«(i') = -1^ (c-'^"(*-*") - 1) e-'^"*°e'^'"aL 



h.c. 



and the double integrals of the commutators yield 



(F.6) 



■^m Sm 



t .t' 
di' / dt" 

to J to 

t ^t' 

X / dt' / 



dt"sin(5™(A'-<")- (/''-^(^')) 



dt' 



cos I - 



•cos((5„(A'-to)- (^''^-^^^'^ 



lit - to) cos 



- sin (S,nit - to) - ((^(■'' - (^(^■))) - sin ((/.^'^ - <^(j')) 

(F-7) 

Note that in the time evolution operator corresponding 
to 7i(t) the terms sin (v?(*) - ip'^^^) = - sin (</j(j) - <^(*)) 
cancel each other. 
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Appendix G: Canonical Transformation 



The commutators for read 



The unitary operator of the canonical transformation 
has the form 



C/c := 

and is apphed to 
n := (eat + ^ 



(G.l) 




(G.2) 

Here, the constants ^ e C and k G M. 

We use the Baker-Campbeh-Hausdorff formula from 
section B to do the transformation. The commutators 
for Hi yield 



[Ca^ + Ca, Aa^ - X*a] = ^X* (-at 



a + aa' 



— e*A (d^a — aa' ^ 

2?r . 



Commutators with higher "nesting levels" trivially van- 
ish. Hence, the complete transformation of Hi reads 

n'l := C/e-HiC/t = (Cat + ^d) + ^1 (G.4) 

K 



[—Ko'a, Adt — A*d] = — k (d^ [d, Ad^ — A*d] 
+ [dt, Adt - A*d] d) 
= —K (Ad^ + A*d) 

= -(Cdt + rd) 



(G.5) 



and 



[-Kdtfl, Adt - A*d] = [- (Cdt + Co) , Adt - A*d] 

(G.6) 

where we used equation (G.3). Higher order terms in the 
expansion trivially vanish again. The complete transfor- 
mation of reads 

n'2 UcW^ = -na^'a - (^d^ +e*d) - (G.7) 

K 

Hence, the transformation of the full Hamiltonian 
reads. 



H' ='H'i+'H'2 = — i - Kd^d. 



(G., 



We want to stress that the canonical transformation is 
exact in this case. 
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